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ABSTRACT 


The  true  submersible,  capable  of  extended  operating  periods  without 
surfacing,  has  long  been  desired.    Present  day  submarines  must  surface 
in  order  to  charge  batteries.    Celestial  navigation,  therefore,  is  available 
as  a  by-product  of  normal  operating  procedures.    With  the  advent  of  the 
submarine  capable  of  extended  periods  of  operation  submerged,  it  will  be 
necessary  to  provide  a  method  of  navigation  that  will  not  require  the  sub- 
marine to  surface,  if  full  advantage  is  to  be  taken  of  her  other  character- 
istics. 

The  Instrumentation  Laboratory,  Massachusetts  Institute  of  Technol- 
ogy, proposes,  as  a  solution  to  the  navigation  problem,  the  application  of 
inertial  techniques  employing  gyros  to  provide  the  inertial  reference  and 
stabilized  platforms  to  provide  position  indication.    The  problem  is  thus 
one  of  measuring  the  change  in  a  direction,  that  of  the  true  vertical,  with 
respect  to  inertial  space. 

Inertial  techniques  are  being  applied  to  air  navigation.      For  these  ap- 
plications, it  is  assumed  that  the  carrying  vehicle  departs  from  a  stable 
base,  where  a  well  defined  vertical  and  azimuth  direction  are  available 
for  "lining-up"  the  system  prior  to  departure.   A  somewhat  different  prob- 
lem exists  for  a  submarine  underway  at  sea.    Here,  there  are  no  well 
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established  references  available;  hence,  the  boundary  conditions  are 
largely  unknown. 

A  study  of  possible  systems  for  applying  inertial  methods  to  sub- 
marine navigation  is  currently  underway  at  the  Instrumentation  Labora- 
tory, Massachusetts  Institute  of  Technology.    A  report  (15)   has  been 
issued  on  the  initial  phase  of  the  study,  in  which  only  undamped  systems 
employing  ideal  components  were  considered.    A  self -settling  system, 
however,  must  employ  damping,  which,  together  with  component  per- 
formance, will  affect  the  accuracy  of  the  system;  hence,  the  complete 
study  will  consider  these  effects  on  system  performance.    Furthermore, 
the  resulting  system  must  have  a  period  approaching  the  84.  4-minute 
Schuler-tuning  condition  if  it  is  to  be  effective.    This  fact  implies  that 
the  transient  effects  due  to  disturbances  imposed  upon  the  system  will 
have  a  duration  measured  in  terms  of  hours.    It  would  therefore  be  de- 
sirable to  have  some  means  available  for  predicting  the  errors  resulting 
from  such  disturbances. 

This  thesis  suggests  that  the  mean-square  error  in  position  can  be 
predicted  for  any  particular  time  during  the  duration  of  the  transient  by 
using  the  power  spectral  density  of  the  disturbance  and  a  function  defined 
as  the  transient  transfer  function,  and  that  this  method  may  be  applied  in 
the  design  stage  to  arrive  at  component  performance  specifications. 

For  illustration,  the  equations  for  one  proposed  inertial  system,  us- 
ing nonideal  components  which  in  themselves  generate  uncertainties  or 
disturbances,  have  been  derived  and  the  application  of  the  proposed  meth- 
od for  evaluating  the  resulting  errors  is  shown. 

It  is  regretted  that  lack  of  time  and  non -availability  of  equipment  pre- 
vented a  more  thorough  study  of  the  problem.    However,  from  the  work 
completed,  it  appears  unsafe  to  design  a  system  on  the  basis  of  steady - 
state  servoanalysis  techniques  where  the  solution  of  the  problem  during 
the  duration  of  the  transient  is  also  of  primary  importance.    This  conclu- 
sion is  based  upon  the  fact  that  the  response  of  the  third-order  system 
studied  to  a  forcing  input  e^^^  made  an  initial  excursion  which  was  up 
to  several  times  as  large  as  the  steady-state  response.  Furthermore, 
the  ratio  of  the  initial  excursion  to  the  steady-state  response  increased 
with  an  increasing  frequency,  whereas  the  time  at  which  the  maximum  ex- 
cursion occurred  did  not  appear  to  be  any  well  behaved  function  of  the  fre- 
quency. 
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OBJECT 

To  investigate  the  possible  application  of  a  self- 
contained  inertial  system  to  the  navigation  of  sub- 
merged vehicles. 
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CHAPTER  1 
DISCUSSION  OF  ESfERTIAL  NAVIGATION 

INTRODUCTION 

The  problem  of  navigation  of  a  submerged  vehicle  is  of  interest 
in  view  of  proposed  submarine  operating  characteristics  that  enable 
extended  submerged  periods.    These  characteristics  will  make  it  nec- 
essary to  navigate  from  the  submerged  vehicle  for  considerable  per- 
iods of  time,  with  some  degree  of  accuracy,  and  without  reference  to 
external  information.    Solution  of  the  problem  of  long-range  all-weath- 
er aircraft  navigation  has  been  attempted  by  inertial  techniques.    An 
attempt  to  apply  similar  techniques  to  the  navigation  of  underwater 
vehicles  is  only  now  being  investigated.  ^     '' 

Six  systems  for  inertial  navigation  have  been  proposed  and  studied 
in  the  undamped  form,  using  ideal  components.  ^     '    This  thesis  pro- 
poses to  investigate  the  effects  on  the  performance  of  one  of  these  sys- 
tems when  nonideal  components,  which  introduce  uncertainties  such 
as  gyro  drift,  integrator  drift,  etc. ,  are  included  in  the  system,  and 
further,  to  extend  the  study  to  the  damped  system. 

INERTIAL  NAVIGATION 

Navigation  is  the  science  of  guidance  of  a  vehicle  over  the  Earth, 
over  or  through  its  waters,  and  through  its  atmosphere.    This  guid- 
ance involves  a  two -fold  function  of  indication  and  control.    Indication 
is  the  comparison  of  the  indicated  position  of  the  vehicle  with  some 
reference  fixed  in  the  Earth.    This  comparison  may  be  used  to  obtain 
the  distance  to  a  destination  and  the  direction  the  vehicle  must  move 
to  reach  that  destination.    It  has  thus  effectively  produced  the  control 
signal. 
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Piloting,  which  is  the  simplest  form  of  guidance,  is  ignored  in 
this  discussion,  because  it  is  assumed  that  no  optical  line  of  sight  is 
available.    The  substitution  of  electro -magnetic  beams  for  optical 
lines  of  sight,  and  the  substitution  of  time  intervals  for  distance  in- 
formation, leave  the  essentials  of  piloting  in  the  guidance  problem. 
This  guidance  system  is  subject  to  interference,  either  natural  or 
man  made.    The  desire  to  minimize  interference  has  stimulated  the 
search  for  a  self-contained  guidance  system  operating  on  such  natural 
earth  phenomena  as  are  effectively  interference  proof. 

Inertial  navigation  can  be  regarded  as  a  precision  dead -reckoning 
method  in  which  measurements  are  made  with  respect  to  inertial  space 
rather  than  with  respect  to  the  water  surrounding  the  vehicle.    Conven- 
tional celestial  navigation  is  a  form  of  partial  navigation  in  which  stars 
represent  the  inertial  space  references;  but  the  requirement  for  sub- 
merged operation  precludes  the  consideration  of  celestial  systems.    Full 
inertial  systems  use  gyroscopes  for  inertial  reference. 

Physical  Factors  Available  for  Guidance 

The  natural  phenomena  associated  with  the  Earth  that  offer  possibi- 
lities for  use  with  a  self-contained  guidance  system  include  data  from  the 
medium  surrounding  its  vehicle,  the  Earth's  electric  field,  the  Earth's 
magnetic  field,  and  the  Earth's  gravity  field.    References  (13)  and  (1) 
discuss  each  of  these  possibilities  in  general  and  conclude  that  the  Earth's 
gravity  field  appears  to  be  the  best-suited  source  of  data  for  a  self- 
contained  guidance  system. 

The  mass  of  the  Earth  has  associated  with  it  a  gravitational  field, 
which  field  is  essentially  radial  to  the  center  of  the  Earth  and  has  no 
poles.    The  direction  of  this  gravitational  field  could  be  readily  deter- 
mined by  a  plumb  bob  if  the  Earth  were  non-rotating.    Due  to  a  basic 
law  of  physics,  namely,  the  principle  of  equivalence  in  the  general 
theory  of  relativity,  that  gravitational  mass  is  equivalent  to  inertial  mass, 
the  daily  rotation  of  the  Earth  produces  a  centrifugal  acceleration  that 
also  affects  the  plumb  bob.    The  vector  resultant  of  the  Earth's  gravita- 
tional field  and  the  centrifugal  force  field  is  defined  as  the  Earth's  gravity 
field,  and  the  direction  of  this  essentially  radial  gravity  field  is  unique 
at  any  given  point  on  the  Earth  (Fig.  1-1). 
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Fig.  1-1   Gravitational  field  intensity  and  gravity  field  intensity  (the  vertical). 


The  Earth  itself  is  subject  to  the  equivalence  of  gravitational  and  in- 
ertial  effects.   The  associated  centrifugal  force  field  causes  the  Earth  to 
bulge  at  the  equator,  giving  it  a  nearly  ellipsoidal  or  spheroidal  shape. 
The  figure  that  a  fluid  body  with  the  mass  distribution  and  daily  rotation 
of  the  Earth  would  have  is  defined  as  the  geoid.   Attempts  are  made  to  rep- 
resent the  geoid  by  an  associated  reference  ellipsoid  for  the  basis  of 
geodetic  measurements  and  map  data.    The  surface  of  the  geoid  is  rep- 
resented by  mean  sea  level,  which  is  everywhere  an  equipotential  sur- 
face of  the  Earth's  gravity  field.   Because  the  geoid  does  not  have  a 
consistently  smooth  surface,  the  vertical  is  not  in  general  parallel  to  the 
normal  to  the  reference  ellipsoid  at  the  same  position.    The  angle  between 
the  vertical  and  the  ellipsoid  normal  is  called  the  deflection  of  the  vertical. 
This  deflection  is  generally  less  than  0.30  milliradian  ;  and  generally,  the 
largest  marine  deflections  occur  in  regions  of  rapid  changes  in  ocean  depth. 

In  general,  variations  in  the  deflection  of  the  vertical  can  be  consid- 
ered to  be  quasi-static  at  submarine  Speeds,  and  can  be  expected  to  have 
no  dynamic  effect  on  the  operation  of  submarine  inertial  navigation  systems. 

Basic  Principles  of  Inertial  Navigation 

The  vertical,  which  is  the  direction  of  the  specific  force  of  gravity, 
is  unique  at  any  point  on  the  Earth  and  is  the  basis  of  the  astronomical 
position  of  that  point.   Astronomical  position  is  a  direction  relative  to  the 


Earth  and  has  two  components: 

1.  Astronomical  latitude,  the  angle  between  the  equatorial  plane  and 
the  vertical  at  the  located  point.    This  quantity  is  unique  because  the 
Earth's  polar  axis  has  a  unique  spatial  direction. 

2.  Astronomical  longitude,  the  angle  between  the  plane  through  the 
Earth's  polar  axis  containing  an  arbitrary  vertical  (e.  g. ,  that  at  Green- 
wich) and  the  plane  through  the  polar  axis  containing  the  located  vertical. 
This  quantity  itself  is  not  unique  because  of  its  arbitrary  reference,  but 
difference  of  longitude  is  unique. 

Astronomical  position  data  will  not  agree,  in  general,  with  geodetic 
position  data  (on  which  maps  are  based)  due  to  the  before  mentioned  de- 
flection of  the  vertical.    However,  position  data  determined  by  vertical 
indication  and  corresponding  map  data  should  agree  in  the  case  of  small 
land  masses  whose  positions  have  been  plotted  from  astronomical  data. 

The  fundamental  problem  associated  with  inertial  navigation  is  then 
the  indication  of  the  vertical  inside  a  self-contained  vehicle  that  is  under- 
going additional  linear  accelerations.    It  is  this  problem  which  has  been 
considered  in  detail  in  Reference  (15). 

Theory  of  Vertical  Indicators 

The  simplest  vertical  indicator  is  the  plumb  bob  or  pendulum.    A  pen- 
dulous element  is  defined  as  a  system  free  to  rotate  about  a  point,  so  that 
it  assumes  a  preferred  orientation  under  the  influence  of  impressed  ex- 
ternal forces.    If  the  suspension -point  is  not  accelerated  with  respect  to 
the  Earth,  the  pendulous  element  can  provide  a  precise  indication  of  the 
true  vertical,  i.  e. ,  the  direction  of  gravity.    If  the  suspension-point  is 
accelerated  with  respect  to  the  Earth,  the  pendulous  element  tends  to  in- 
dicate the  apparent  vertical,  i.  e. ,  the  direction  of  the  resultant  specific 
force. 

When  the  pivot  of  the  pendulous  element  is  accelerated,  a  reaction 
torque  acts  on  the  center  of  the  mass  of  the  element,  producing  a  torque 
which  causes  the  element  to  lag  behind  the  pivot.    To  study  the  effects  of 
an  acceleration  without  introducing  the  effects  of  the  mass  of  the  element, 
it  is  convenient  to  discuss  the  motion  in  terms  of  force  per  unit  mass,  i.  e. , 
in  terms  of  the  inertia  reaction  specific  force  associated  with  the  accelera- 
tion.   An  inertia  reaction  specific  force  has  the  magnitude  of  its  associated 
acceleration  (though  expressed  in  different  units)  and  is  opposite  in  direc- 


tion.    Figure  1-2  is  a  force  diagram  for  an  accelerated  pendulous  ele- 
ment in  a  gravitational  field,  indicating  that  the  direction  of  the  result- 
ant specific  force,  i.  e. ,  the  apparent  vertical,  is  the  vector  resultant 
of  the  gravitational  specific  force  and  the  inertia  reaction  specific 
force. 
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Fig.  1-2    Specific   forces  acting  on  a  pendulous  element. 


To  obtain  the  specific  forces  acting  on  a  pendulous  element  during 
a  navigational  problem,  it  is  necessary  to  know  the  associated  accelera- 
tions.   The  linear  acceleration  of  a  point  with  respect  to  a  reference 
coordinate  system  may  comprise  five  components  with  respect  to  a 
coordinate  system  that  is  moving  relative  to  the  reference  coordinate 
system.    However,  as  far  as  marine  navigation  is  concerned,  only  two 
components  are  significant. 

1.  The  tangential  (horizontal)  acceleration  associated  with  the  angu- 
lar acceleration  of  the  true  vertical  as  the  vehicle  accelerates  over  the 
Earth's  surface. 

2.  The  centripetal  acceleration  parallel  to  the  projection  of  the 
Earth-radius  on  the  equatorial  plane,  due  to  the  Earth's  daily  rotation; 
this,  added  to  the  gravitational  field  intensity,  gives  the  gravity  field  in- 
tensity. 

When  a  vehicle  moves  over  the  nearly  spherical  earth,  the  vertical 
associated  with  its  instantaneous  position  rotates  approximately  geocentri- 


cally  with  respect  to  the  Earth.    On  a  ship  moving  at  constant  speed  along 
a  great  circle  over  the  Earth,  a  pendulum  would  have  to  rotate  with  re- 
spect to  the  Earth  at  an  effectively  constant  angular  velocity  in  order  to 
indicate  the  vertical  continuously.    Once  the  transient  stages  were  over, 
no  torque  would  be  required  to  keep  the  pendulum  along  the  vertical  dur- 
ing this  rotation.    K  the  vehicle  were  to  accelerate,  by  changing  either 
speed  or  heading,  a  torque  would  be  required  to  maintain  the  pendulum  on 
the  vertical,  which  in  turn  would  be  geocentrically  accelerating. 

Schuler  ^     '  pointed  out  in  1923  that  accurate  indication  of  the  vertical 
during  periods  of  acceleration  could  be  realized  by  suitable  tuning  of  the 
natural  period  of  the  pendulum.    Whenever  the  pivot  of  a  pendulum  is  ac- 
celerated, the  center  of  mass  of  the  pendulum  tends  to  lag  behind  the  piv- 
ot with  respect  to  inertial  space.    At  the  same  time,  this  acceleration 
causes  the  true  vertical,  associated  with  the  pivot,  to  accelerate  geocen- 
trically with  respect  to  inertial  space.    These  considerations  lead  to  the 
following  observation: 

If  a  pendulum  initially  hangs  vertically,  it  will  remain  along  the  ver- 
tical if  its  angular  acceleration  about  its  pivot  equals  the  geocentric  ang- 
ular acceleration  of  the  vertical.  (The  accelerations  are  with  respect  to 
inertial  space.) 

These  two  accelerations  become  equal  for  a  distributed -mass  pendu- 
lum when  the  ratio  of  the  pivot -center  of  mass  separation  to  the  square 
of  the  radius  of  gyration  of  the  pendulum  equals  the  reciprocal  of  the  radius 
of  the  Earth.    When  operation  is  in  the  Earth's  gravity  field,  and  the  pen- 
dulum is  undamped,  this  condition  gives  a  pendulum  with  a  natural  period 
of  approximately  84  minutes.    For  a  shorter -period  pendulum  the  rotation- 
producing  torque  on  the  pendulous  element  is  stronger,  so  the  pendulum 
lags  the  true  vertical.    For  a  longer-period  pendulum  the  torque  is  weaker 
and  the  pendulum  leads  the  true  vertical. 

An  84-minute  period  is  very  long  compared  to  the  periods  of  physical 
pendulums  ordinarily  encountered.    In  fact,  it  is  improbable  that  a  simple 
(concentrated-mass)  pendulum  or  a  distributed -mass  pendulum  could  be 
constructed  with  this  period,  although  the  84-minute  gyropendulum  just 
fails  of  being  realizable.    However,  in  the  practical  case,  a  condition  anal- 
ogous to  the  Schuler-tuning  of  a  pendulum  can  be  realized  in  a  particular 
kind  of  closed-loop  system  in  which  a  short-period  pendulum  tracks  the 


the  specific  forces.    A  platform,  controlled  on  the  basis  of  the  pendulum 
output  data,  then  becomes  an  equivalent  Schuler  pendulum. 

From  the  foregoing,  it  is  deduced  that  a  vertical  indicating  system  is 
conceivable  which,  while  not  a  pendulum,  utilizes  the  Schuler -tuning  idea, 
so  that  a  torque  caused  by  the  vehicle's  acceleration  acting  on  the  system 
will  be  absorbed  by  the  inertia  reaction  torque  required  to  keep  a  refer- 
ence line  in  the  system  on  the  vertical. 

Reference  (15),  from  which  the  foregoing  section  was  obtained,  has 
elaborated  upon  this  theory  and  set  forth  six  representative  inertial  sys- 
tems for  position  indication. 

Fundamentals  of  Position  Indication 

Astronomical  position  measurement  involves  the  determination  of  the 
essentially  geocentric  angles  between  indicated  and  reference  verticals. 
It  is  proposed  to  determine  these  angles  by  using  a  stabilized  platform  as 
an  essential  feature.    The  stabilization  loop  receives  specific  force  data 
from  accelerometers  or  pendulums,  and  delivers  these  data  to  gyro  units 
and  servo  drives,  which  in  turn  drive  the  platform  toward  a  horizontal 
position.    The  specific  force  data  are  modified  and  filtered  to  control  the 
system  dynamics.    Part  of  this  modification  consists  of  Schuler -tuning 
the  stabilization  loop.    Two  basic  methods  of  utilizing  the  indicated  ver- 
tical are  then  available  to  obtain  the  vehicle's  essentially  geocentric  angle 
of  travel.    One,  called  here  analytical  integration,  involves  the  integra- 
tion of  accelerometer  outputs  or  of  gyro  inputs,  with  the  accelerometers 
and  gyro  units  mounted  on  the  stable  platform.    The  other  method,  called 
here  geometric  integration,    involves  the  direct  comparison  of  the  indicated 
vertical  with  a  reference  vertical.    Both  methods  have  been  outlined^  ^^ 
in  terms  of  six  representative  systems  utilizing  ideal  components  and  un- 
damped systems,  namely: 

(1)  Open-Chain  Integration  of  the  Angular  Velocity  of  the  Indicated 

Vertical. 

(2)  Open-Chain  Double  Integration  of  the  Angular  Acceleration  of 

the  Indicated  Vertical. 

(3)  Direct  Double  Integration  of  Acceleration. 

(4)  Electromagnetic  Composition  of  Tangential  and  Vertical  Com- 

ponents of  Celestial  Longitude  Rate. 


(5)  Geometric  Integration  Using  Angular  Velocity  of  the  Indicated 

Vertical  with  Respect  to  Inertial  Space. 

(6)  Geometric  Integration  Using  Pre -aligned  Inertial  Reference. 
The  first  system  wUl  now  considered  in  some  detail  using  nonideal 

components  in  both  the  damped  and  undamped  systems,  with  special  em- 
phasis on  the  problem  of  determining  a  possible  method  for  setting  up 
performance  specifications  for  components.  It  should  be  noted  here  that 
the  normal  methods  of  frequency  analysis  do  not  apply  directly  where  the 
system  does  not  reach  steady  state  during  its  normal  operating  cycle,  or 
where  the  transient  period  of  the  system  is  of  primary  interest.  Such  is 
the  case  in  the  operation  of  an  inertial  position  indication  system. 

Open-Chain  Integration  of  the  Angular  Velocity  of  the  Indicated  Vertical 

Figure  1-3  is  a  functional  diagram  of  a  single -degree -of -freedom 
system  for  indicating  the  vertical.    Figure  1-4  is  a  functional  diagram  of 
the  essential  operating  components  of  a  latitude  indicator  which  uses  an 
open-chain  single  integration  of  the  angular  velocity  of  the  indicated  ver- 
tical to  arrive  at  the  indicated  latitude.    Figure  1-4  corresponds  to  Fig. 
B-1  in  Appendix  B.    Note  that  Fig.  B-1  is  primarily  the  basic  vertical  in- 
dicator of  Figs.  A -3  and  1-3,  with  the  angular  velocity  integration  tapped 
off  the  vertical  indicator  loop  at  the  angular  velocity  signal  that  calls  for 
processing  the  gyro  unit,  so  that  the  position  integrator  output  is  the  in- 
dicated change  in  latitude.    The  vertical  indicator  is  undisturbed  by  the 
presence  of  this  integrator,  as  far  as  the  loop  dynamics  are  concerned. 
Hence,  the  same  Schuler -tuning  conditions  applies  to  this  system  as  ap- 
plied to  the  basic  vertical  indicator  of  Fig.  A -3. 

A  two -degree -of -freedom  system  constructed  along  these  lines  would 
require  an  x-axis  accelerometer  in  addition  to  the  y-axis  accelerometer 
shown  in  Fig.  1-4,  as  well  as  x-axis  components  corresponding  to  the 
orientational  control -signal  generating  system,  gyro  unit,  and  controlled 
member  drive.    In  addition,  the  x-axis  system  requires  that  its  angular 
velocity  signal  be  modified  by  secant  of  latitude,  in  order  to  give  longitude 
rate  data  to  its  position  integrator.    The  position  integrator  in  the  x-axis 
system  would  indicate  longitude  difference  with  respect  to  some  arbitrary 
reference.    The  maintenance  of  the  y-and  x-axes  along  a  parallel  of  latitude 
and  a  meridian  respectively  must  be  handled  by  an  auxiliary  orientation 
system,  e.g.,  a  gyrocompass.  (15) 
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The  gyro  units  ideally  maintain  their  spin  axes  fixed  with  respect  to 
inertial  space.    To  maintain  the  controlled  member  on  which  they  are 
mounted  tangential  to  the  Earth  in  its  daily  rotation  therefore,  the  x-axis 
gyro  must  be  precessed  by  a  correction  torque  that  keeps  the  platform  up 
with  Earth  rate.    The  calculations  of  Appendix  B  consider  this  torque  to 
be  generated  by  a  special  precession  current  applied  to  the  torque  gener- 
ator of  the  gyro  unit. 

In  Appendix  B,  it  is  shown  that,  if  the  controlled  member  has  the 
equation  of  motion  of  a  Schuler-tuned  pendulum,  the  correction  to  the  in- 
dicated position  (that  is,  the  instantaneous  difference  between  the  true  and 
indicated  positions)  has  a  similar  equation  of  motion,  except  for  the  ap- 
pearance of  certain  initial  conditions  as  additive  constants  of  integration 
in  the  indicated  position.    Specifically,  the  correction  to  the  indicated 
position  at  any  time  is  the  correction  to  the  indicated  position  at  the 
start  of  the  problem,  plus  the  present  correction  to  the  indicated  vertical. 
The  indicated  position  is  thus  Schuler -periodic  with  respect  to  boundary- 
matching  terms  associated  with  the  vertical  indication.    The  above  descrip- 
tion holds  directly  with  ideal  components  throughout  the  system.    Position 
indication  will  be  modified  somewhat  however,  where  uncertainties  are 
present  in  the  system,  as  is  shown  in  Appendices  A  and  B. 
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CHAPTER  2 
PROPOSED  METHOD  FOR  ARRIVING  AT  COMPONENT  SPECIFICATIONS 
UNCERTAINTIES  IN  COMPONENTS 


QUANTITY 

SMOOTHED 
QUANTITY 

J^^--z^rr^^^''^\^^^ 

QUANTITY 

q       -s^"" 

^^^^""Z^Z- — — \" 

> 

1 

QUANTITY 
INACCURACY 

ROUGHNESS 

REFERENCE 
QUANTITY 

CORRECT 
QUANTITY 

(Rghns)q 

a 

(^^1^^ 

(E)q,s., 

(Rgh)q  =  q  -  q    =  (U)q 

?>? 

^-^ 

SMOOTHED 

(R)qJ         '   ' 

QUANTITY 

■^' 

QUANTITY 

ERROR 

(E)q(s„„=a    -0    , 

^(corr) 

CORRECT 
LEVEL 

^(sm)    ^(corr) 

(i)q  =q  -  q,^„„, 

TIME  t  *- 

Fig,  2-1  Graphical  definitions  of  inaccuracy,  roughness,  and  error,  for  a 
quantity  varying  with  time. 

Figure  2-1  is  a  graphical  definition,  by  Draper  ^  ',  of  inaccuracy, 
roughness,  and  error  for  a  quantity  varying  with  time.    The  relationships 
shown  are  summarized  below  for  convenience. 


{R)q 


=  {U)q     =  q  - 


^(sm) 


(E)q 
(I)q 


(sm)  ,  ^(sm)  ~  ^(corr) 


^  "  ^(corr) 
=  (U)q  -  (E)q 


also 


(sm) 


(C)q 


(sm) 


=  -(E)q 


(sm) 


(2-1) 

(2-2) 

(2-3) 
(2-4) 

(2-5) 


In  the  remainder  of  this  thesis,  when  the  error  and  correction  in  q 
are  referred  to,  it  will  be  taken  to  mean  (E)q/      \  and  (C)q(gjjj). 
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The  ideal  components  used  in  the  inertial  systems  described  in  Re- 
port R-9  ^     '  have  been  replaced,  for  the  purpose  of  this  study,  with  com- 
ponents which  are  no  longer  ideal  in  that  the  individual  components  now 
are  considered  as  introducing  spurious  signals  into  the  system  by  their 
operation. 

For  the  purposes  of  this  study,  uncertainty  in  the  output  of  an 
integrator  has  been  treated  as  an  equivalent  spurious  input  to  the  inte- 
grator ^     '.   The  specifications  for  an  integrator  in  which  the  process 
has  no  specified  maximum  value  may  be  written: 

<U)l(out)(int)     S     ^'=^/l(i„)(spec)   "t-  (2-6) 

That  is,  the  specification  is  such  as  to  require  that  the  uncertainty  of 
the  output  quantity  of  the  integration  process  be  less  than  a  percentage 
of  the  magnitude  of  the  integral  of  some  specified  value  of  the  input 
quantity  to  the  integrating  system,  over  the  interval  of  integration.   In 
this  connection,  it  must  be  assumed  that  the  uncertainty  of  the  integration 
process  will  not  be  truly  random  over  the  period  of  integration  because 
of  the  fact  that  the  integrator  may  be  subject  to  inaccuracies,  such  as 
drift  and  frictional  losses,  which  do  not  lend  themselves  to  correction; 
therefore,  the  average  uncertainty  of  the  integrating  process  will  be 
different  from  zero  over  the  period  of  integration.    Similarly,  the  un- 
certainty in  the  output  of  an  amplifier  may  not  be  a  truly  random  pro- 
cess over  a  given  interval  because  of  amplifier  drift,  and  the  average 
uncertainty  will  therefore  be  different  from  zero  over  the  interval. 
From  the  foregoing,  it  appears  that  specifications  similar  to  those 
placed  upon  the  integrator    may  be  placed  upon  the  amplifier,  i.e., 
that  the  uncertainty  in  the  output  shall  be  less  than  some  percentage  of 
the  magnitude  of  the  output  due  to  a  specified  value  of  the  input  quantity 
taken  over  a  specified  period  of  time. 

The  foregoing  may  be  represented  mathematically  as  in  Equation 
(2-7)  where 

m^n  <        ^(^"^P)^n>    gout)    ^  .,  „. 

^U)q(out)(amp)      -  J^^  g(in)(spec)'  ^^"^^ 

or  its  equivalent 

(U)q  <         g(put)(spec)        over  the  specified  time 

(out)  (amp)  N  interval. 
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It  should  be  noted  that,  in  both  the  integrator  and  the  amplifier,  the 
specification  implies  that  the  input  signal  has  undergone  an  ideal  inte- 
gration or  amplification,  and  that  superimposed  upon  the  output  is  an  in- 
accuracy which  is  due  wholly  to  an  uncertainty;  specifically,  that  the  pro- 
cesses of  integration  and  amplification  are  carried  through  without  error. 
Uncertainties  in  gyro  units  and  controlled  member  drives,  pendulum  units, 
and  accelerometer  units  as  developed  by  the  Instrumentation  Laboratory 
are  generally  susceptible  to  separation  into  two  categpries:   those  due  to 
uncertainties  in  electronic  signals  applied  as  inputs  to  the  torque  summing 
member  or  derived  as  outputs  from  the  position  of  the  torque  summing 
member;  and  those  due  to  torque  uncertainties  acting  on  the  torque  sum- 
ming member,  such  as  mechanical  unbalances,  bearing  friction,  and 
undesirable  reaction  torques  applied  to  the  torque  summing  member  by 
the  microsyn  signal  units,  etc.    The  subject  has  been  treated  extensively 
in  several  references,  (10),  (7),  (6),  and  (16),  and  is  of  such  scope  that 
no  attempt  will  be  made  to  analyze  it  here,  beyond  the  general  observa- 
tion that  the  results  of  open -chain  tests  would  seem  to  indicate  that  the 
powerful  methods  of  statistical  analysis  might  be  necessary  in  order  to 
determine  the  suitability  of  units  of    various  designs  for  inertial  system 
applications.    This  method,  which  is  essentially  that  recommended  in  later 
sections, again  would  treat  the  units  as  ideal,  with  uncertainties  superim- 
posed on  the  output  of  the  ideal  unit;  hence,  preliminary  specifications  may 
be  laid  down  following  principles  to  be  developed  later. 

In  the  previous  section,  uncertainties  are  defined,  and  the  manner  in 
which  uncertainties  may  arise  in  various  components  are  discussed.    Now, 
it  is  necessary  to  consider  the  effect  of  uncertainties  when  the  components 
from  which  they  arise  are  employed  in  closed-loop  systems.    In  particular, 
performance  of  components  in  a  closed-loop  system,  as  for  instance  the 
vertical  indicating  loop  of  Fig.  l-3,affects  the  dynamics  of  the  loop,  there- 
by determining  the  precision  of  indicating  the  vertical.    Components  which 
operate  in  open-chain  fashion  outside  the  loop,  as  for  instance  the  position 
indicating  integrator  of  Figs.B-1  and  B-2,  have  direct  additional  effects  on 
the  precision  of  the  system.     Hence,  the  open-chain  specification  method  , 
as  proposed  in  previous  paragraphs,  as  well  as  the  method  for  arriving  at 
performance  specifications  for  closed-loop  components,  is  of  direct  ap- 
plication. 
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By  the  assumption  that  q/4n\  is  operated  on  by  a  particular  component 
ideally,  but  that  the  component  itself  then  contributes  an  uncertainty  in  the 
output  which  is  superimposed  upon  the  true  output,  we  have  in  effect  sepa- 
rated the  total  output  ihto  the  sum  of  two  outputs,  one  of  which  is  the  true 
output,  and  the  other,  the  portion  of  the  total  output  due  to  the  uncertainty. 
If,  then,  it  were  possible  to  predict  by  some  simple  means  the  magnitude 
of  the  uncertainty  and  the  way  in  which  it  varies  with  respect  to  some  vari- 
able, say  time,  then  it  should  also  be  possible  to  determine  what  effect  the 
uncertainty  would  have  upon  the  performance  of  the  system. 

Unfortunately,  the  very  term  uncertainty  implies  the  nature  of  the  dif- 
ficulty. We  cannot  by  any  simple  means  determine  the  necessary  informa- 
tion on  the  behavior  of  the  uncertainties.  However,  by  statistical  study  of 
the  behavior  of  components  under  various  operating  conditions,  it  may  be 
possible  to  determine  the  probable  form  of  the  uncertainties.  In  this  re- 
port, it  is  considered  that  the  results  of  such  a  study  would  be  presented  in 
the  form  of  power  spectral  densities,  which  would  then  be  used  to  enter  the 
response  curves  for  a  given  system,  in  order  to  determine  the  effect  on  sys- 
tem performance  of  a  component  with  an  output  uncertainty  of  a  certain  spec  - 
tral  density. 

Since,  in  a  position  indication  system,  we  are  dealing  with  a  dynamical 
system  of  second -order  or  higher,  and,  since  our  systems  are  lightly  damped, 
the  system  may  have  normal  modes  or  resonant  frequencies  that  can  be  ex- 
cited by  forcing  functions  of  the  same  frequency.   In  particular,  Wrigley  has 
shown  ^     '  that  an  undamped  vertical  indication  system  of  the  type  described 
later  in  this  thesis  is  insensitive  to  a  forcing  functionwith an  84.4-minute  period. 

The  mathematical  expression  of  the  weighting  function  describes  these 

k  D  t 
modes,  which  are  found  to  be  the  terms  of  the  form  t  e^i   into  which  the 

weighting  fimction  can  be  broken.^  ^In  a  complicated  system,  it  may  be  very 
difficult  to  obtain  an  explicit  mathematical  statement  for  W(t),  and,  in  de- 
sign problems,  the  amount  of  work  involved  in  optimizing  system  parameters 
for  the  presence  within  the  system  of  certain  exciting  frequencies  may  be  ex- 
cessive. A  further  discussion   of  the  weighting  function  is  included  in  later 
sections.   On  the  other  hand,  as  shown  in  Appendix  D,  it  may  be  a  relatively 
simple  matter,  using  machine  methods,  to  determine  the  steady -state  trans- 
fer function,  Y(jw),  the  weighting  function  and  the  transient  transfer  func- 
tion, F(w,t),  which  is  defined  and  discussed  later;  and  from  these  to  con- 
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struct  charts  from  which  system  specifications  can  be  prepared. 

The  purpose  of  the  following  sections  is  to  present  in  detail  the 
manner  in  which  such  information  may  be  utilized  either  in  preparing 
specifications  or  in  predicting  the  performance  of  the  system. 

Again,  it  should  be  emphasized  that  the  necessity  for  evaluating 
position  information  at  all  times  the  system  is  in  operation  requires 
a  detailed  consideration  of  the  transient  which  persists  for  a  consid- 
erable length  of  time  because  of  the  long  natural  period  of  the  system, 
and  the  long  system  characteristic  time. 

THE  TRANSIENT  TRANSFER  FUNCTION 

Consider  some  arbitrary  system  performance  function  that  may 
be  represented 

VUL   .  (P-'-q)  .  (3-8) 

^(in)  P    ^agP^+ajp+a^ 

The  transfer  function  Y(p)  is  frequently  defined  as  the  expression  on 
the  right  with  initial  conditions  zero,  which  is  also  the  Laplace  trans- 
form of  the  weighting  function,  where  the  weighting  function  is  defined 
as  the  response  of  the  system  to  a  unit  impulse  with  the  system  initially 
at  rest.  ^^^ 

In  the  usual  frequency  response  studies  of  a  dynamical  system,  the 
steady-state  response  to  a  periodic  function  is  the  item  of  interest.    In 
an  inertial  navigation  system,  however,  because  of  the  long  system 
periods  involved,  the  system  may  and  probably  will  be  placed  in  use  be- 
fore the  steady-state  has  been  achieved.    Furthermore,the  effects  of  sys- 
tem disturbances  are  present  as  transients  for  such  a  length  of  time  that 
a  study  of  the  transient  response  of  the  system  to  periodic  forcing  functions 
appears  necessary. 

In  order  to  study  this  transient  response,  we  will  define  a  new  para- 
meter such  that,  if  a  system  represented  by  a  performance  function  of 
the  form  of  Equation  (2-8)  is  subjected  to  a  forcing  function  of  the  form 
cos  wt  +  j  sin  wt  =  e^^  ,  the  response  may  be  represented  as  follows: 

jwt 


e-* 


^      '  system 


F(w,t)eJ^ 


Fig.  2-2.  System  performance  function. 
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where,  for  the  steady-state,  Lim  F(w,  t)—»-Y(w),  the  transfer  function 

t  — ^co 
for  the  system  in  which  Y(w)  will  now  be  defined  as  the  steady -state 

response  of  the  system  to  the  forcing  function  e^^f  .    The  F(w,  t)  then 
represents  the  transfer  function  of  the  system  during  the  transient  por- 
tion of  the  response,  and  is  the  parameter  in  which  we  are  interested; 
moreover,  it  may  be  obtained  by  subjecting  the  system,  initially  at 
rest,  to  forcing  inputs  of  the  form  e^      for  a  series  of  frequencies  and 
obtaining  the  time  response  of  the  system.    |f(w,  t)|    thus  obtained  may 
be  plotted  vs.  time  for  various  frequencies.    From  these  plots,  cross 
plots  of  |f(w,  t)|      for  various  times  vs.  frequency  can  also  be  prepared 
(Fig.  3-17).    The  manner  of  obtaining |F(w,t) I      is  of  interest  and  is 
covered  in  Appendix  (D). 

Now,  postulating  that  any  disturbance  to  which  the  system  is  likely 


to  be  subjected  can  be  assigned  a  power  spectral  density  G(w),we  may 

2, 


write  for  the  steady-state:  ^  ' 


G(out)(^)  =|y(w)|  'G(i„)(w).  (2-9) 

Analogously,  for  the  transient  period, 

G(out)("'')    =|F(w,t)|2G(j„)(w).  (2-10) 

Further,  it  can  be  shown  that  the  ensemble  mean  square  value  of  the  out- 
put quantity  at  any  time  t^  for  any  arbitrary  input  with  power  spectral 
density  given  by  G/^  v(w),  is  the  integral  of  the  product  of  the  spectral 
density  of  the  input  and  the  square  of  the  transient  transfer  function  at 
time  t^  between  the  limits  w  =  0  to  w  =  oo,  that  is, 

MS  output  at  time{t^)=f  |F(w,tj)|  ^G/.  x(w)  dw.  (2-11) 

It  should  be  pointed  out  here  that  the  ensemble  mean  square  may  be 
defined  as  the  mean  of  the  squared  end -point  values  of  q(t)  at  time  tj,  for 
a  number  of  trials;  whereas  the  time  mean  square  commonly  used  in  com- 
munications applications  is  the  time  averaged  mean  square  of  q(t)  over  a 
specified  time  interval.    Since  we  are  interested  in  the  output  of  our  sys- 
tem at  any  specified  time  rather  than  the  average  output  over  some  arbi- 
trary time  interval,  it  follows  that  the  ensemble  mean  square  of  the  out- 
put is  the  quantity  of  interest. 

Now  supposing  |f(w,  t)|      to  be  determined  by  some  such  means  as  that 
described  in  Appendix  (D);  a  plot  might  then  be  constructed  presenting 
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|f(w,  t)  I     vs.  frequency,  as  shown  in  Fig.  3-18.    It  will  be  noted  that 
by  using  this  plot,  and  assuming  that  G/jj^x(w)  may  be  represented  as  a 
constant,  the  integration  of  Fig.  3-17  could  be  performed  as  the  product 
of  G/j  x(w)  and  the  area  under  the  curve  of  |  F(w,  t)  |    .    Or  more  general- 
ly, given  plots  of  |  F(w,  t)  |     and  spectral  density  of  a  forcing  function 

G/._x(w)  vs.  w  for  any  system  under  consideration,  neither  curve  being 

^    '  I  I  2 

necessarily  representable  analytically,  where  |  F(w,  t)  |     is  plotted  for 

time  tj  and  spectral  density  is  time  invarient,  the  integration  of  Equation 

(2-11)  may  be  performed  by  any  of  several  standard  methods  of  graphical 

integration.  ^     ^ 

An  examination  of  Fig.  3-18  reveals  that  the  open-chain  uncertainty 
specifications  previously  discussed  are  insufficient  to  specify  component 
performance  in  a  closed-loop  system,  since  no  account  is  taken  of  pos- 
sible frequency  components  present  in  the  uncertainties.    Fig.  3-18,  for 
example,  shows  that  frequency  components  in  the  vicinity  of  FR  =  1  will 
have  far  greater  effect  than, for  instance,  the  same  magnitude  of  uncer- 
tainty with  frequency  components  in  the  vicinity  of  FR  =  4.    However,  by 
making  certain  assumptions,  it  is  possible  to  use  the  information  obtained 
from  a  study  of  |  F(w,  t)  |     to  arrive  at  prelihiinary  component  specifications. 

If  the  uncertainty  output  of  a  component  is  assumed  to  be  a  purely  ran- 
dom process,  i.  e. ,  the  uncertainty  has  the  characteristic  of  "white  noise", 
the  power  spectral  density  of  the  uncertainty  is  then  constant  over  an  in- 
finite band  of  frequencies  from  zero  to  infinity.    However,  I  F(w,  t)  I     may 
approach  zero  if  the  system  with  which  it  is  associated  has  an  upper  cut- 
off frequency,  that  is,input  frequencies  beyond  a  certain  upper  value  of 
frequency  will  have  no  discernible  effect  on  the  system  output.    All  the 
systems  considered  in  this  report  have  this  upper  frequency  cut-off  char- 
acteristic, i.  e. ,  the  systems  behave  as  low -pass  filters,  and  the  assump- 
tion of  "white  noise"  may  therefore  be  of  use. 

It  should  be  noted,  however,  that  an  assumption  of  "white  noise"  of  a 
given  amplitude,  is  one  of  the  worst  which  can  be  imposed  upon  a  system, 
since  the  integration  then  involves  all  frequencies  with  equal  weight.    It 
seems  reasonable  to  assume  that  the  actual  operating  component  will  give 
rise  to  uncertainties  whose  frequencies  will  be  predominantly  centered 
within  certain  frequency  bands,  with  negligible  components  outside  these 
bands.    If  then  the  maximum  spectral  density  in  the  predominant  bands 
were  to  be  applied  to  the  system  as  "white  noise"  of  the  same  amplitude, 
the  resulting  system  performance  would  be  worse  than  the  actual  by  the 
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value  of  the  integral  outside  the  predominant  frequency  bands. 

The  above  reasoning  suggests  that  preliminary  specifications  might 
be  written  on  the  basis  of  the  open-chain  performance  specifications  if 
it  is  assumed  that  the  uncertainties  thus  obtained  are  the  worst  possible 
"white  noise"  cases.    Then,  any  improvement  which  might  occur  be- 
cause of  more  favorable  frequency  distributions  may  be  considered  as 
a  welcome  gratuitous  improvement  in  system  performance. 

System  specifications  for  position  indicating  systems  are  generally 
written  as  requiring  that  indicated  position  be  within  specified  limits  of 
true  position  during  specified  time  limits.    The  procedures  of  the  fore- 
going yield  information  on  the  mean -square  errors  in  system  indication 
due  to  uncertainties  in  individual  components.   It  can  be  shown  that  the 
sum  of  the  magnitudes  of  these  individual  mean -square  errors  is  equal 
to  the  mean-square  error  in  the  system  output  when  due  regard  is  given 
to  the  possibility  that  the  individual  terms  may  be  complex.    Hence,  the 
procedure  followed  in  the  remainder  of  this  report  is  to  analyze  the  sys- 
tems on  the  basis  that  each  component  in  turn  is  nonideal,  the  remaining 
components  in  the  system  being  considered  perfect.    The  information 
thus  obtained  should  be  useful  in  evaluating  which  components  contribute 
the  most  to  system  error,  and  in  determining  lower  marginal  require- 
ments which  must  be  met  by  components  to  be  considered  for  use  in  the 
system. 

THE  WEIGHTING  FUNCTION 

The  weighting  function  (W(t))  may  be  defined  as  the  inverse  trans- 
form of  the  term  on  the  right  of  Equation  (2-8).    For  further  discussion 

of  this  subject,  the  reader  is  referred  to  any  good  text  on  servomechan- 

^9\  (2) 
isms.    ^  '  ^  '    The  weighting  function  for  any  system  can  be  determined 

by  subjecting  the  system  to  a  unit  impulse,  i.e., 
where 


*<"  =  (PF),y,,,„  (Uj)t, 


(u.  )t  =  unit  impulse. 

Conversely,  since  the  unit  impulse  is  the  derivative  of  the  unit  step, 
the  weighting  function  is  also  the  derivative  of  the  response  of  the  system 
to  a  unit  step  input. 

The  weighting  function  has  the  characteristic  that  it  relates  the  pre- 
vious inputs  to  the  system  to  the  present  output  by  the  integral 
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^(out) 

(t) 

Y_<.(tn) 

(t- 

-T) 

W(7-) 

dr, 

where 

r  = 

t  - 

-h 

• 

(2-12) 


Hence,  a  knowledge  of  the  weighting  function  enables  one  to  form  an  esti- 
mate of  the  transient  errors  present  in  the  output  of  the  system  following 
a  forcing  function  which  is  applied  over  a  given  period  of  time  and  then 
goes  to  zero,  such  as  the  transient  errors  in  an  inertial  system  caused 
by  maneuvers  of  the  ship.    This  conclusion  follows,  if  it  is  assumed  that 
any  given  maneuver,  occurring  over  a  period  of  time,  produces  a  forced 
error  in  the  output  of  the  systein.    If  the  ship  then  returns  to  a  steady 
course  and  speed  at  time  (tj),  the  forcing  function  is  no  longer  being  ap- 
plied to  the  system,  but  the  error  caused  by  the  maneuver  remains  in 
the  system  at  time  (t)  as  a  transient  related  to  the  forced  error  at  time 
(tj^)  by  the  value  of  the  integral  in  Equation  (2-12).    The  above  example 
is  a  specialized  application  of  the  general  theory  by  which  a  series  of 
impulse  or  step  functions  may  be  used  to  approximate  any  arbitrary  in- 
put to  the  system,  and  thus  the  normal  response  of  the  system  to  the  in- 
put may  be  determined  from  the  relationship  of  Equation  (2-12). 
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CHAPTER  3 
ANALYSIS  OF  VERTICAL  INDICATORS 

GENERAL 

The  differential  performance  equation  for  the  position  indicating  sys- 
tem being  considered  in  this  thesis  is  derived  in  Appendix  B.    The  differ- 
ential performance  equations  for  the  associated  vertical  indicators,  both 
damped  and  undamped,  are  derived  in  Appendix  A.    Equation  (B-7)  shows 
that  the  correction  to  the  indicated  vertical, [(C)V]/^  ^x,  enters  directly  as 
one  of  the  terms  in  the  correction  to  the  indicated  position,  [(C)P]/.   -y 
Consequently,  in  the  future,  when  references  are  made  directly  to  the 
performance  of  the  vertical  indicator,  the  reader  may  infer  that  refer- 
ence to  position  indication  is  implied  as  well. 

In  this  chapter,  we  will  consider  first  the  undamped  vertical  indicator, 
and  finally  the  damped  vertical  indicator.    The  concept  of  a  transient  trans- 
fer function  as  introduced  in  Chapter  2  will  be  used  as  an  aid  in  evaluating 
the  effects  of  certain  uncertainties  upon  the  vertical  indication.    Since  the 
results  of  any  uncertainty  spectral -density  analyses  on  the  various  com- 
ponents have  not  been  available  to  the  authors,  if  indeed  such  studies  ex- 
ist in  usable  form,  only  qualitative  estimates  will  be  made  as  to  the  rela- 
tive performances  required  of  the  various  components. 

UNDAMPED  VERTICAL  INDICATOR 

Figure  3-1  presents  in  summary  form  the  equations  of  |F(w,t)  |     for  the 
undamped  vertical  indicator.    Figure  3-2  is  a  semi-log  plot  of  Equations 
(3-5)  and  (3-6)  vs.  frequency  ratio  for  a  time  of  one -half  hour.    Photo- 
graphs of  transient  transfer  functions  of  damped  vertical  indicators  are 
shown  in  Fig.  3-3.    These  damped  systems  are  to  be  discussed  later  in 
this  chapter.    It  will  be  recalled  that  the  mean-square  errors  obtained  in 
the  proposed  method  are  the  ensemble  mean  squares,  and  consequently  rep- 
resent the  errors  at  a  particular  instant  of  time.     Further,  by  the  previous 
definitions,  the  time  interval  is  measured  from  the  instant  the  forcing 
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airen  the  parforaaact  aqtatton  of  a  tytteii  for  which  the  responae  to  a  forcing  input  of  the  font  e^  f  la  dealred,  and  reeallln«  that  e'  f* 
eoa  Vft  *  J  aln  W.t,  the  reaponae  of  the  aysten  to  e  f  la  given  by  the  eomplez  avin  of  the  aeparate  reaponaea  of  the  ayaten  to  eoa  W.t  and  ain  W^t. 
We  can  apply  the  foregoing  to  the  undaaped  vertical  Indicator  with  the  differential  equation  Qq.  A-32)^irtileh  uy  be  written 


[P»  *  WJJ  [(0)»](,_,)  -  S(,^^_^^^.o,)  <l(i„)o  *  ^q(i„)i.O,)  "   "^dn)!.  {i-32) 

where  q/,_)T  and  q/,^),  correspond  to  the  appropriate  quantltiea  appearing  in  the  lero 

and  firat  order  derlvatlvea  reapectively  of  equation  (A-32)  is  i^pendiz  A. 
Denote  by  '/.w..  \(t)i  the  reaponae  of  the  ayatea  to  q.,  w  s  •  eoa  W.t  where  the  aubseript  (n)  la  the  order  of  the  derivatlTe  aaaociated  with 
*(ln)(n)'  "^  *'  '(n)(eln)^*''  ****  "•?">"••  <>'  *''•  ayatem  to  <l/jj,w  \   "  aln  W^t.  Then  taking  initial  eonditlona  equal  to  aero,  the  aolutlon  for  Bq. 
(A-32)  to  the  forcing  function  e  f  may  be  written  aa  the  eum  of  the  four  teraa  below. 

eoa  W,t  -  eoa  W._t 
^o)(coa)(*)-    wJ-W,»   "  0-« 

J'(o)(aln)(»>  ■  i    (  Wj-W,0  ('^  V  -  7^    'l"  "nl*)  0"*' 

'(l)(coa)(»)  '    {  Vrii^  )  Cf  •!»  V  -  *n«  -^  W  <«' 

•"a)(.in)(*'  =  4w7^)  (-"•  V*  "•*».»'  o-*> 

Becalllng  the  definition  of    |r{w,t)|.  Chapter  2,  and  making  appropriate  aubstitutlone  following  the  method  of  Appendix  D,  We  obUln 

l'(v.t)l»^^  -    l^(o)(co.) (*)!'*'•'  '(o)(ai«)<»)l'  •  (wX^-i^)    '     [<"•  V  -  "•  "nX*)'  *  ^"^  V  "  ^  •^»  W> 

(3-5) 

IT(w.t)l"(,)  .    I'd)  (eoa)  (t)!"*    1^  '(l)(ain)*l'  "  (i^-^Ff^"  ^^"^  "°  "'^  "  """  "°  ''"'*^'  *  ^°"  ""**  "  "'  ''''^'^'       ^^^ 

er,  the  tlae  aolntlona  for  the  aero  order t era  of  the  danped  Tertleal  indicator  of  Equation  (A-38)  nay  be  written 

1/2 


'(o)(coa)(»>-  ^^^,  .(i,^)'3   [(ww^  .  Aj'  .  ,^j  a  .  p«),,]  •     *|[(5|^r*w,»][(w„J-w,»)»..w„JW»w,»]|   -^^^V 


^^ : w,-i ^s TT -rr    e"^  '  na*  aln  (W^Jl  -  (DB)»  ♦  A,)t 

W^  Vl  -  (IB)"   [[(c|-  -  (DH)W^J  ♦  WJ  (1  -  (DH)»)]  [[W_^J  (2(DH)»-  1)  +  W,»]»*  I*  W^J  {DR)»  (1  -  (DH)»)]J 

"  (3-7) 

'^"^•"  '*'  ■■  [v*(ct)l(-^;-4)-*w„j(i-(i.)f  ""  W-.']KJ-v)'-Va(-)'vj)     '"'"^'  * "^^ 

,  I1/2  -(DH)W  t        ^ 

r5 ^^.-i^ j\         e      "  aln  (»„  Jl  -  (DB)"*  A  -A^)t . 

r-  (DR)W  )  ♦  W  ;  (1  -  (DH)»)J[[W  J  (2(DH)a-l)  ♦  W,a]»  *  H  M*    (DB)»  (1  -  (DB)a)jJ 


0-8) 

2(DB)W.  W, 


"na   "f 


,   _,  -1^1  -  (DB)»   ,  .1   ''n.V^  -  t"")'    ,  -1  "nS  (PH)Vl  -  (DB)» 
^2  =*"   -TSh^-*"      1  .(j^),    -*"    WJ(2(DB)».1)*W,4 


A3  .  tan-lL^H  . 
The  corroaponding  |T(w,t)  |,  ,  for  Bquationa  (3-7)  and  (3-8)  la  not  included  here,  but  the  reader  can  readily  verify  that  the  reaulting  function 
aooawhat  conplez,  and  that  r(w,t)L j  la  not  a  particularly  well  behaved  function  of  the  forcing  frequency. 

Fig.  3-1   Summary  of  time  solution  for  vertical  indicators.     . 
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Fig.  3  -2   Transient  transfer  function  squared  of  undamped  vertical  indicator 
for  time  t  =  1/2  hour. 
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functioii  e^^f   is  applied  to  the  instant  in  time  under  consideration. 
Also,  although  strictly  speaking  the  frequency  form  of  the  transfer  func- 
tion is  not  defined  for  an  unstable  system  ^  %  the  frequency  form  of  the 
transient  transfer  function  does  exist  by  the  definition  in  Chapter  2. 

The  area  under  |  F(w,t)|  /q^  plotted  on  linear  coordinates  would  rep- 
resent the  integral    f   |F(w,t)|  ^G/^j^xCw)   dw  with  G/^^^x  =  0,  0<w<  oo,  and 
is  less  than  the  area  under  the  curve  for    F(w,t)j/,s  plotted  in  the  same 
fashion.    Hence,  for  "white  noise"  for  which  G/.   x(w)  is  constant  for  all 
frequencies,  the  term  containing  the  first  derivative  produces  the  greater 
error  in  vertical  indication. 

The  first  derivative  term  is  associated  with  the  uncertainty  in  gyro 
drift  and  with  the  inaccuracy  in  the  speed  compensation  current.   Con- 
versely, the  uncertainties  in  the  accelerometer  unit  and  the  integrator 
within  the  vertical  indication  loop  are  associated  with  the  zero -order 
term.    Therefore,  the  performance  of  the  gyro  unit  and  the  measurement 
of  water  speed  will  probably  be  more  critical  than  the  performance  of  the 
accelerometer  unit  and  the  integrator.    This  conclusion  will  be  especially 
true  if  the  spectral  densities  of  the  imcertainties  have  a  wide  band  width 
such  that  the  dips  in  |  F(w,t)  |     for  the  first-order  term  do  not  have  an  ap- 
preciable effect  in  reducing  the  associated  errors. 

The  one -half -hour  time  was  chosen  as  representative  of  possible  sit- 
uations in  which  a  damped  system  would  not  be  used.    These  situations 
would  be  those  in  which  the  period  of  operation  would  be  short  with  re  - 
spect  to  the  time  required  to  damp  out  transient  disturbances,  and  hence 
no  advantage  would  be  gained  from  the  greater  complexity  of  the  damped 
system. 

DAMPED  VERTICAL  INDICATOR 

Next,  we  consider  the  damped  vertical  and  position  indication  system. 
The  parameters  used  are  those  that  provide  the  optimum  position  indica- 
tion with  respect  to  the  angular  velocity  of  the  true  vertical.   The  result- 
ing parameters  are  summarized  in  Table  C  -I. 

In  Reference  (8),  the  optimization  was  done  on  the  opposite  basis  from 
that  used  in  Appendix  C;  specifically,  the  acceleration  error  was  eliminated 
and  the  jerk -error  minimized.   It  is  interesting  to  note  that  the  attenuation 
of  errors  due  to  accelerations  under  the  conditions  of  Reference  (8)  pro- 
ceeds at  6  decibels  per  octave  with  increasing  frequency,  whereas 
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by  using  the  present  methods,  acceleration  errors  are  attenuated  at  18 
decibels  per  octave. 

Figures  3-4  through  3-15  present,  respectively,  the  steady-state 
sinusoidal  responses  of  the  vertical  indicator,  the  step  responses,  and 
the  impulse  responses  or  weighting  functions,  for  DR  =  0. 1,  0.3,  and 
0.  5.    As  defined  in  Chapter  2,  |  F(w,  t)|is  an  analogous  function  to|Y(w)|, 
the  transfer  function;  and  hence  the  area  under  the  curve  of  |y(w)|      has 
the  same  utility  for  the  steady-state  as  the  area  under  the|F(w,  t)|      has 
for  the  transient  period.    Hence,  from  a  study  of  the  steady-state  re- 
sponse curves  and  a  comparison  with  Equation  (A -38),  we  may  conclude, 
by  reasoning  similar  to  that  applied  in  the  case  of  the  undamped  system, 
that  the  uncertainties  associated  with  the  gyro  unit  and  water-speed  com- 
pensation will  probably  be  most  critical.    Proceeding  in  a  similar  man- 
ner, we  may  arrive  at  the  following  listing,  which  is  essentially  a  list- 
ing of  components  in  direct  order  of  the  excellence  of  performance  re- 
quired, but  without  regard  to  the  sensitivities  which  may  be  involved. 

(1)  (a)    Gyro  unit 

(b)   Water -speed  compensation 

(2)  Accelerometer  unit 

(3)  (a)    Direct -channel  amplifier 

(b)  Earth-rate  compensation 

(c)  Feedback  amplifier 

(d)  High-frequency  integrator 

(4)  Low-frequency  integrator 

It  should  be  noted  that,  in  the  damped  vertical  indication  system, six 
sensitivities  appear.    These  six  sensitivities  are  related  in  Derivation  Sum- 
mary A-U  by  only  four  equations.    Thus,  we  may  impose  two  more  condi- 
tions on  the  sensitivities, such  that  we  arrive  at  six  independent  equations 
in  the  six  sensitivities.    It  is  apparent,  then, that  we  may  be  afforded    some 
leeway  in  adjusting  sensitivities  such  as  minimizing  the  errors  due  to  some 
particularly  objectionable  component  in  the  system. 

No  mention  has  been  made  in  this  chapter  of  the  effect  on  position  in- 
dication of  uncertainties  in  the  open-chain  integrator.    However,  as  shown 
in  Appendix  B,  the  mean-square  error  in  position  due  to  the  open-chain  in- 
tegration may  be  developed  in  a  manner  analogous  to  that  applied  above. 

The  remainder  of  this  chapter  is  devoted  to  a  discussion  of  the  methods 
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Fig.  3-4   Steady-state  sinusoidal  response  of  three  subsystems  of  damped 
vertical  indicator,  DR  =  0.1. 
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Fig.  3-5   Steady-state  sinusoidal  response  of  three  subsystems  of  damped 
vertical  indicator,  DR  =  0.3. 
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Fig.  3-6   Steady-state  sinusoidal  response  of  three  subsystems  of  damped 
vertical  indicator,  DR  =  0.5. 
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Fig.  3-7    Steady-state  sinusoidal  response  of  two  subsystems  of  damped 
vertical  indicator,  DR  =  0.1. 
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Fig.  3-8    Steady-state  sinusoidal  response  of  two  subsystems  of  damped 
vertical  indicator,  DR  =  0.3. 
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Fig.  3  -9    Steady  state  sinusoidal  response  of  two  subsystems  of  damped 
vertical  indicator,  DR  =  0.5. 
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Fig.  3-10  Step  response  of  damped  vertical  indicator,  DR  =  0.1. 
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Fig.  3-11    Impulse  response  of  damped  vertical  indicator  (derivative  of 
step  response),  DR  =  0.1. 
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Fig.  3-12   Step  response  of  damped  vertical  indicator,  DR  =  0.3. 
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Fig.  3-13   Impulse  response  of  damped  vertical  indicator,  (derivative  of 
step  response),  DR  =  0.3. 
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Fig.  3-14   Step  response  of  damped  vertical  indicator,  DR  =  0.5. 


37 


[P]  =  p3  +    a2p2  +  QiP   +  Qc 


STEP   INPUT 


[PF]  =  p       Mag.  Ix 


[PF]=-pp    Mag.lx 


[PFJ  =  -^       Mag.  Ix 


[PF]  =  p-      Mag.  Ix 


[PF]  =  -*^^'       Mag.lx 


Fig.  3-15  Impulse  response  of  damped  vertical  indicator,  (derivative  of 
step  response),  DR  =  0.5. 
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used  and  the  results  obtained  in  studying  the  transient  transfer  func- 
tion for  the  damped  vertical  indicator. 

For  the  study  of  the  damped  vertical  indicator,  it  was  attempted 
to  apply  the  method  described  in  Appendix  D,  using  a  high-speed,  gen- 
eral purpose  analog  computer. 

Two  identical  systems  (Fig.  3-16a)  were  placed  on  the  simulator 
and  the  forcing  functions  cos  w^t  and  sin  w^t  were  applied,  one  to  each 
system.    The  output  of  one  system  was  then  taken  to  the  x-axis  of  an 
oscilloscope.    The  resulting.  Lissajous  pattern  represented  a  plot  of 
the  sine  response  of  the  system  vs.  the  cosine  response,  from  which 
it  was  proposed  to  obtain  the  transient  transfer  function,  |  F(w,  t)  |  . 

It  was  intended  that  the  sine  and  cosine  functions  be  developed  in 
an  harmonic  oscillator  (Fig.  3-16b)  as  a  part  of  the  set-up.    The  oscil- 
lator would  then  oscillate  at  the  frequency  determined  by  the  w,  setting 
when  excited  by  a  step  input  from  the  computer  master  generator. 

Theoretically,  such  an  oscillator  would  be  possible  using  an  un- 
damped second-order  system,  and  the  sine  and  cosine  outputs  plotted 
against  each  other  would  produce  a  circle  on  the  oscilloscope.    However, 
due  to  nonlinearities  or  other  imperfections  in  the  computer,  it  was 
found  that  the  amplitude  of  the  oscillations  of  both  terms  tended  to  change 
in  magnitude.    Consequently,  it  was  found  necessary  to  introduce  either 
positive  or  negative  damping  into  the  system,  depending  upon  the  fre- 
quency, to  maintain  a  constant  amplitude. 

The  damping,  in  turn,  affected  the  sine  and  cosine  outputs  differ- 
ently; hence,  it  was  necessary  to  include  the  "scope-attenuator"  in  each 
system  to  adjust  the  resulting  ellipse  to  a  circle.    The  use  of  an  addi- 
tional attenuator  in  this  position  will  not  introduce  distortions  in  the  sys- 
tem responses,  since  all  signals  are  attenuated  the  same  amount  for  a 
given  setting. 

The  "input  gain  control"  was  put  into  the  input  to  the  harmonic  os- 
cillator to  provide  a  calibrated  gain  control  in  order  to  increase  or  de- 
crease the  system  responses  to  usable  size.  In  tests,  changes  in  set- 
ting of  this  control  seemed  to  produce  no  distortion  in  the  circular  dis- 
play of  the  sine  vs.  cosine  pattern,  and  hence  it  may  be  assumed  that 
the  amplitudes  of  both  were  changed  equally. 
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Fig.  3-16   Block  diagrams  of  high-speed  analog      computer. 
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The  remainder  of  the  set-up  is  ciearly  apparent,  the  remaining  gains 
and  attenuators  being  applied  merely  to  provide  more  reasonable  potentio- 
meter settings. 

In  calibrating  the  oscilloscopes,  one  of  which  was  used  as  a  monitor- 
scope,  and  the  other  as  a  photo-scope,  the  necessary  damping  ratio  was 
first  set  into  the  harmonic  oscillator.    Then,  with  K/.  \  and  K    set  to  con- 
venient values,  the  oscilloscope  amplifiers  were  used  to  set  in  a  two-inch 
diameter  circle  on  the  monitor -scope.    Then,  with  K/^^^x  reduced,  and  all 
other  settings  maintained  constant,  a  two-inch  circle  was  adjusted  on  the 
photo-scope  using  the  oscilloscope  amplifiers.    The  shift  in  K/.  \  in  the 
calibrations  was  made  so  as  to  take  advantage  of  the  linear  amplifiers  on 
the  oscilloscope  to  amplify  small  pictures,  rather  than  alternatively  ap- 
plying large  forcing  signals  with  the  attendant  danger  of  distortion  due  to 
over -loading  the  computer. 

After  calibration,  the  oscilloscope  amplifiers  were  left  undisturbed. 
For  each  frequency  thereafter,  a  two-inch  circular  plot  of  the  sine  vs.  the 
cosine  was  obtained  on  the  monitor -scope  by  manipulation  of  the  damping 
ratio  in  the  oscillator  and  the  scope  attenuators,  K  .    The  input  gain  was 
set  at  a  convenient  value  for  the  above  operation  at  each  frequency  and 
recorded,  as  was  the  gain  used  for  photographing  each  response,  in  order 
that  the  responses  might  be  reduced  to  a  common  reference.    Calibration 
was  completed  by  recording  the  two-inch,  or  unit-reference  circle  photo- 
graphically, in  order  to  provide  a  direct  reference  for  reducing  the  data. 
Using  the  above  methods,  photographs  of  the  Lissajous  pattern  responses 
of  each  of  the  sub-systems  were  made  at  various  frequencies.    It  was  ap- 
parent, however,  that  the  highest  forcing  frequencies  which  might  be  ap- 
plied corresponded  to  a  frequency  ratio  of  about  five,  because  of  the  large 
attenuation  in  the  system  at  higher  frequencies  and  the  consequent  limita- 
ations  imposed  by  the  noise  level  in  the  computer.    Thus  the  data  could  not 
be  extended  to  the  significant  frequency  region  for  which  FR  >  16.    The  data 
were  then  reduced  from  the  Lissajous  patterns,  of  which  those  shown  in  Fig. 
3-3  are  representative.    The  time  dots  on  these  patterns  signify  4 -hour 
intervals.    An  examination  proved  these  data  to  be  unreliable,  since  steady- 
state  amplitudes  thus  obtained  failed  to  agree  with  steady-state  amplitudes 
computed  by  other  means.    As  a  consequence,  no  effort  will  be  made  to 
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present  any  quantitative  results. 

Alternatively,  a  qualitative  summary  of  the  results  of  a  study  of  the 
Lissajous  patterns  from  sub-system  to  sub-system  and  damping  ratio 
to  damping  ratio  is  presented  below.    The  points  in  the  summary  occur 
conListently  and  are  considered  significant  in  the  transient  study. 

(a)  An  excursion  of  the  response  occurs  that  is  greater  than  the 
steady-state  response.    In  addition,  the  initial  excursion  may  or  may  not 
be  oscillatory. 

(b)  The  response  then  returns  toward  the  steady-state  in  an  oscil- 
latory manner,  the  frequency  of  the  oscillations  apparently  increasing 
with  increasing  frequency. 

(c)  The  ratio  of  the  initial  excursion  to  the  steady-state  response 
apparently  increases  with  an  increase  in  frequency. 

(d)  The  time  in  hours  at  which  the  maximum  excursion  occurs  is 
apparently  a  function  of  both  damping  ratio  and  frequency  for  any  parti- 
cular system,  and, for  the  systems  studied,  varied  from  six  to  forty 
hours. 

(e)  The  time  of  reaching  the  maximum  excursion  does  not  appear  to 
be  an  orderly  function  of  frequency  for  a  particular  damping  ratio  and  sys- 
tem.   That  is,  for  increasing  frequencies,  the  time  at  which  the  maximum 
excursion  occurs  may  first  show  an  increase  and  then  a  decrease.    From 
the  data  obtained,  no  reliable  conclusion  seems  possible  as  to  the  effect 

of  damping  ratio  on  the  time  for  reaching  the  maximum  excursion. 

(f)  The  nature  of  the  responses  seems  to  indicate  the  inadvisability 
of  designing  a  system  purely  on  the  basis  of  conventional  steady -state 
analysis  techniques,  if  the  response  during  the  transient  period  is  of  in- 
terest. 

The  previous  statements  are  not  considered  by  the  authors  to  be  a 
condemnation  of  the  method.    It  is  suggested  that  additional  work,  which 
could  not  be  undertaken  by  the  authors  because  of  a  combination  of  lack 
of  time  and  non -availability  of  equipment,  might  yield  interesting  results. 
Conservatively,  these  results  would  show  low-frequency  behaviors  of  the 
sub-systems  for  the  vertical  indicator.    It  would  then  remain  to  determine 
the  relationship,  if  any,  that  exists  between  the  behavior  at  high  and  low 
frequencies. 

A  method  proposed  by  Dr.  J.  H.  Laning,  Jr. ,  of  the  Instrumentation 
Laboratory,  Massachusetts  Institute  of  Technology,  would  approximate 
the  high-frequency  transient  transfer  function  using  the  response  of  the 
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system  to  impulses.    This  method  will  be  investigated  at  some  la*^er  date. 

Another  method,  suggested  by  Dr.  Sidney  Lees,  also  of  the  Instru- 
mentation Laboratory,  would  make  use  of  a  special  applicatic    of  the  roots - 
locus  method  of  servoanalysis  to  reduce  the  amount  of  labor  involved  in 
calculating  a  point  by  point  response  of  the  transient. 

Still  another  method,  for  which  equipment  was  not  available,  would 
again  use  the  analog  computer.    Here,  the  transient  transfer  function 
would  be  presented  as  a  function  of  time  by  squaring  and  adding  the  re- 
sponses of  the  two  systems,  taking  the  square  root  of  the  sum,  and  ap- 
plying the  resulting  signal  to  one  axis  of  an  oscilloscope,  the  other  of 
which  is  the  time  axis. 

The  remaining  material  in  this  chapter  is  presented  not  with  the  idea 
of  showing  any  quantitative  information,  but  by  using  the  data  as  obtained 
to  show  the  manner  in  which  the  theory  might  be  applied. 

Figure  3-17  is  a  linear  plot  of  |  F(w,  t)  |      obtained  by  the  methods  of 
Appendix  D.    It  is  from  this  plot  that  it  is  proposed  that  mean-square  er- 
rors be  determined  by  the  integral  of  Equation  (2-11). 

M  S    outputL^^  =      y'^|F(w,tj)|  ^G(in)(w)    dw.  (2-11) 

1  Q  I 

Figure  3-18  is  the  more  familiar  semi-log  plot  of  Fig.  3-17,  on  which 
is  also  plotted,  but  not  to  scale,  the  steady-state  transfer  function  for  com- 
parison.   Thus,  it  is  seen  that  the  apparent  spikes  in  Fig.  3-17  are  just 
such  as  would  occur  if  the  transfer  function  squared  were  itself  plotted  to 
a  linear  scale.    The  shape  of  the  curve  at  four  hours  is  due  to  the  fact 
that  the  maximum  initial  excursion  has  not  been  reached,  and  further, 
the  particular  system  under  consideration  is  one  for  which  the  time  to 
reach  the  initial  peak  decreases  with  increasing  frequency  over  the  range 
considered.    It  is  interesting  to   note  that  the  remaining  curves  for  12,  24, 
and  36  hours  generally  follow  the  shape  of  the  transfer  function. 

Figures  3-19  through  3-21  are  the  plots  of  the  areas  under  the  curves 
of  Fig.  3-17  for  the  five  sub-systems  defined  by  the  terms  on  the  right  of 
the  simplified  form  of  Equation  (A -38),  shown  below. 
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Fig.  3-21   Ensemble  mean  square  errors  in  damped  vertical  indicator, 

DR  =  0.5. 


48 


t(^)nt,i)=S(q(i„),  C^)   ^   q(in)/Sn)'  ^v)    [FfV)/ 


Pi  ^  s,  -^ 

l(in)'^V^    [P]   ^(in)2         ^^(in)'    ^V^       [P] 


2  ^^^"^    aT 

^  Vn^>Cv)    ^q^n^.   ■'^(qnn^,    C^)    "T^T"  ^(in)(l,0)-^ 


(p^+  a2P) 

+  S.  p    .     ^  ,  (A -38) 

^V)'^V^  [P]  V)(2,l) 


where  [P]  =  characteristic  differential  equation  for 

the  vertical  indicator. 
Following  the  definition  of  Equation  (2-11),  the  plots  of  the  areas 
may  be  assigned  the  special  significance  of  the  mean -square  error  in 
vertical  indication  at  each  time,  t  =  4,  12,  24,  and  36  hours,  for  an 
input  with  power  spectral  density  such  that  G/.  v(w)  =  1  for  1S(FR)<4.  75, 
and  is  everywhere  else  equal  to  zero. 

The  reader  will  note  the  similarity  between  the  shapes  of  the  error 
curves  for  corresponding  systems  at  different  damping  ratios,  and  fur- 
ther, the  apparent  trend  in  the  relative  location  of  the  curves  at  each 
damping  ratio.    However,  considering  the  accuracy  of  the  data  involved, 
it  is  doubtful  that  more  than  a  tentative  conclusion  can  be  drawn  as  to  the 
relative  importance  of  the  sub-systems  as  error  sources  until  further 
work  is  done  along  this  line. 
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APPENDIX  A 
VERTICAL  INDICATION 

BASIC  PRINCIPLES 

The  following  principles  of  vertical  indication  were  set  forth  in  Ref- 
erence(15),  and  are  repeated  here  for  completeness.    Figure  A-1  is  a 
sketch  of  a  three -axis  vertical  indicator. 

The  vector  angle  representing  the  correction  to  the  indicated  verti- 
cal is  expressed  mathematically  as 


[(C)V](t  i)  =  1^^  X  1^^  ,  (A-1) 

where  the  quantities  are  defined  as  in  Fig.  A -2. 

Indication  with  Respect  to  Earth  Axes  —  Accelerometer  Tracking 

The  rate  of  change  of  this  angle  with  respect  to  indicated  position 
axes,  i.  e. ,  in  a  set  of  coordinates  that  consists  of  the  indicated  verti- 
cal and  two  other  mutually  orthogonal  axes  arbitrarily  oriented  in  the 
indicated  horizontal  plane,  is 


[(C)V](t,i)p.=W(gy)t    -W(Ey)j    -W(j,p)Vj   x[(C)V](t^i),  (A .2) 


where  [(C)V]/|^  ^\p      =  rate  of  change  of  indicated  vertical  cor- 

^        rection  with  respect  to  indicated  position 
axes, 

W/gyxj.  =  angular  velocity  of  true  vertical  with  re- 

spect to  the  Earthy 

W/gyv.  =   angular  velocity  of  indicated  vertical  with 

respect  to  the  Earth, 

W/gpvy  =    angular  velocity  of  indicated  position  axes 


with  respect  to  the  Earth  about  the  indicated 
vertical . 
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VECTOR  FORM 

of](C)V]„  i) 
[(C)V](.i)' 


RESULTANT 
SPECIFIC  FORCE 


'(EV)t 

ANGULAR  VELOCITY 

OF  TRUE  VERTICAL 

WITH  RESPECT 

TO  EARTH 


APPARENT 
VERTICAL 


UNIT  VECTOR  ALONG 
INDICATED  VERTICAL 


UNIT  VECTOR  ALONG 
TRUE  VERTICAL 


INDICATED 
VERTICAL 


(t,l) 

CORRECTION  TO  INDICATED 

VERTICAL  WITH  RESPECT  TO 

TRUE  VERTICAL.    (Angle  through 

which  Indicated  Vertical  must 

be  rotated  to  align  it  with 

True  Vertical.) 


AS  A  VECTOR  ANGLE 


[(C)v](t,i)  =   Iv.  X  hp^^-  I^.  X  1        SINCE  [(C)v](t,i)  IS  SMALL 
I       I    t    Sin  i-^x  i 


Fig.  A-2   Geometrical  relationships  associated  with  the  vector  angle  form  of  the  correction  to 
the  indicated  vertical  with  respect  to  the  true  vertical. 
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The  magnitude  of  the  indicated  horizontal  component  of  the  resultant 
specific  forces,  i.  e. ,  the  quantity  sensed  by  an  accelerometer,  is  ex- 
pressed mathematically  in  Equation  (A  -3) . 

(^^)(res)H.  =  I  (^^)(res)H.  H  V^  ^  (^^)(res)| '  (A -3) 

where 

(sf)/  ggXTj   =  indicated  horizontal  component  of  resultant  speci- 

^  fie  force/ 

=  denotes  magnitude  of  the  enclosed  vector  quantity. 

Let  the  system  operate  such  that  the  angular  velocity  of  the  indicated 
vertical  with  respect  to  the  Earth  is  a  function  of  the  magnitude  of  the  in- 
dicated horizontal  component  of  the  resultant  specific  force  and  directed 
so  as  to  rotate  the  indicated  vertical  toward  the  resultant  specific  force 
vector,  namely 


W(EV)i=(P^)(Vi)(f,W)    \\  ^'(^^^(res) 


(A -4) 


where 

(PF)/yjwj  ^x  =  vertical-indicating  performance  function  for  a 
specific  force  input  and  an  angular  velocity  out- 
put. 
Substitution  of  Equation  (A -4)  into  Equation  (A -2),  plus  the  relation- 
ships for  the  resultant  specific  force  from  Derivation  1,  gives 


[(C)V](t,i)p.=W(Ev)t    '^^^\vi){f,Vn    lv.xig-(WEp)E><«E 

-  Wep  X  (W^p  X  Rg)  -  2WjE  X  (Wgp  X  Rg)}    -  W(Ep)y^  X  [(C)V](j^  jj . 

(A -5) 
The  specific  force  components  in  Equation  (A -5)  can  be  conveniently 
resolved  into  components  parallel  to  the  true  vertical  and  perpendicular 
to  the  true  vertical  (the  arbitrary  nature  of  axes  in  the  horizontal  plane 
requires  no  further  subdivision  into  components),  namely, 

^E  =  ~  Vj  ^E  ' 
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^EP  -  V^  ^(EP)V^  -^  W(EV)t  • 

Substitution  of  these  relationships  into  Equation  (A  -5)  gives  the  specific 
force  in  terms  of  components  along  the  true  vertical,  along  the  angular 
velocity  vector  of  the  true  vertical  (which  must  be  perpendicular  to  the 
true  vertical)  and  perpendicular  to  both  of  these,  namely, 

[■(C)^t,i)Pi  =  W(EV)t  -  (PF)(Vl)(f,W)lv,  ><  {\  [e  -  "eW^eV)! 

-  2Be(Wie)  •  W(H;v)t)]  *  W(EV)t«E[W(EP)Vt  *  ^<'^m  ■  Vt'l 

*KWep)  e  ^  ivj]«E )-W(EP)Vj ""  [(c)^](t,  1)  •  (*-6) 

Carrying  out  the  vector  operation  indicated  in  Equation  (A -6),  re- 
grouping terms  and  omitting  the  components  parallel  to  the  indicated 
vertical  as  not  being  applicable  to  the  indication  of  the  vertical    gives 


[<C)V](t,i)p.  ^  (PF)(vi)(f,w)  («f)(res)Vj  [(C)V](t,i)  =  W^^yjt 
"<^*'>(Vi)(f,W)''E<*(EV)t>Pj  "  <*'^>(V1)(£,W)^e[^(EP)Vj 


+  2(Wjg.ly^)  ]  (ly^  X  W(Eyjt)  -  W(Ep)Vj  X  [(C) V](t_  J) ,  (A  .V) 

where 
(««(res) Vt  =  g  -  «EW^EV)t  "  ^Re^Wie"  W(Ev)t)  -  <*  -8) 

=  component  of  resultant  specific  force  along  the  true 
vertical . 


Let  (PF)/yjW£  ^  be  the  following  operator: 

(pr')(vi)(f, w)  =  hvmML  <  "" '  (*-^' 

f 

where     ^  P.  denotes  integration  of  a  function  of  the  indicated  axes  with 
respect  to  the  time.    Substitution  of  Equation  (A -9)  into  Equation  (A -7) 
gives  Equation  (A -10): 
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[(C)V](t,i)p.  +  S(vi)(f,w)      f^^  («^)(res)V,[<C)V](t  i)dt  =  W(E^), 


Pi 


^  2(WjE.  ly/KV.  ^  W(Ev)t)  dt  -  W(Ep)v.  X  [(C)V](t,i).  (A 


10) 


Differentiating  (A -10)  with  respect  to  the  time  gives  the  differential  per- 
formance equation  in  terms  of  functions  of  the  indicated  axes: 


[<C)Vl(t,l)P;  *  hvmM  («'>(res)Vt  [(C)V](t,i) 

=  (1  -S(vi)((,W)  I'E*(^<EV)t>Pj  -  ^(vl)(t,W)  ''e  [^(EP)Vj 

A. 

.  2(WiE.  \)](\^  X  W(Ev)t)  -  (W(Ep)Vj  X  [(C)V](t_  i)p^) ,  (A-11) 

where  it  is  assumed  that  variations  in  R„  are  slow  enough  to  be  treated 

as  quasi-static  and  variations  in  (sf)/_ ^_v,,  are  either  slow  enough  to  be 

ires;  v^ 

quasi-static  or  fast  enough  to  be  ineffective. 
For  Schuler -tuning  let 

S(Vl)(f,W)  =  —  ' 
"s 

where  R^  is  a  set  value  made  as  nearly  equal  as  possible  to  the  magnitude 
of  Rg.    If  it  is  then  assumed  that  R^  =  Rg,  Equation  (A -11)  becomes 

_i^  (^^)(res)V,       _  _  _         _ 

[(^^^3(t,  i)Pi  ^     —^ t^^) V](t,  i)  =  -  [W(EP) v^  -  2  (W(jE) .  1^^)](1^, 

^  W(Ev)t) "  ^^{mv,  ^  t^^)^l(t,  i))  Pi .  (^-12) 

The  terms  W^^p^^^  (i^_  ><  W(EV)t)    ^^^  2  (Wj^.  1^^)(I^^  x  w^^^^^)  are 

(^^)(res)V. 
small  in  comparison  with  .  ,  and  can  be  neglected.    This  can  be 

shown  as  follows:  s 

W/gpvy  ,  the  magnitude  of  the  rate  of  change  of  the  meridian  about 

the  true  vertical  with  respect  to  the  Earth,  is  given  by 
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W/^„x„   =       "     sin  (Co)tan  (Lat) ,  (A -13) 


(EP)Vt"   T^ 


where 


Vg      =  Velocity  of  base , 

Co      =  course  of  base, 

Lat       =  latitude, 

^R 

^EV)t            p        ' 

^(EP)V^  ^ V.  ^  ^(EV)t^ 

-W(EP)V^  W(EV)t, 

(A -14) 


(A -15) 


and 


2(WjE  •  \  (1 V.    ><  W(E v)t)  I  =  2W(E v)t  WiE  ^i^  (Lat) . 


Assume  that  the  base  is  traveling  due  east  at  10  knots  at  70"  latitude: 

=  00.02  min.  of  arc^  (A-16) 


^(EP)V^  ^(EV)t  ^s^(EP)V^  ^(EV)t 


(^^^(res)V^  (^^^(res)V^ 


2WjE  sin  (Lat)  W^^^^^  _    2R^  W^^  sin  (Lat)  W^^^^^ 

^^^^(res)V^  ^^^^(res)V^  ^         arc 


00.  25  min.of 
re. 
(A -17) 


Neglecting  the  above  terms,  Equation  (A -12)  reduces  to 
(sf). 


'(res)V^ 


t^^^^l(t,i)Pi   ^   -^ ^[(C)V](     )  =  -(W(Ep)v.x[(C)V](t^i)). 

^  ^s  ^  (A -18) 

The  term  (W/p,pxy     X  [(C) V]/.   .J  represents  an  orientation  transfer,  not 

a  magnitude  effect,  and  represents  a  secondary  coupling  between  the  x  and 
and  y  systems.    Since  this  will  only  produce  a  "spiraling  in"  effect  as  the 
system  approaches  its  null  position,  it  will  be  ignored  in  subsequent  dis- 
cussions. 
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Equation  (A -18)  now  reduces  to 

-:-^  ^^^^(res)V,     

[(c)v](t,i)p  - -- — -myxt,i)    =0-  (A-19) 

"s 

INSTRUMENTATION  OF  A  SINGLE -AXIS  OF  THE  VERTICAL  INDICATOR 
OF  FIGURE  A-l(UNDAMPED) 

Definitions  of  Symbols  Used  in  the  Following  Derivations 

a^  Coriolis  acceleration 

a^g  Acceleration  of  base  with  respect  to  Earth 

ayg  Acceleration  of  base  with  respect  to  inertial  space 

p.  Indicated  position  (latitude  or  longitude)  of  vehicle 

(Pj)q  Indicated  position  of  vehicle  at  start  of  problem 

P.  True  position  of  base 

[(C)P]/j.  ^x  Position  correction 

[(C)P]/.   .X  Position  correction. at  start  of  problem 

[(C)V]/j.  jx  Correction  to  the  indicated  vertical 

[(C)V]/x  i\  Correction  to  indicated  vertical  at  start  of  problem 

[(C)V]/x.     ^  Correction  to  reference  vertical 

[(C)V]/j  j.\q  Correction  to  reference  vertical  at  start  of  problem 

6/    X  Output  signal  voltage  of  device  indicated  by  subscript 

G  Gravitational  field  vector 

g  Gravity  field  vector 

H  Angular  momentum  of  gyro  unit  rotor 

i/  \  Output  current  of  device  indicated  by  subscript 

M  Torque  on  gyro  gimbal 

R  =  Rg  Radius  of  the  Earth,  a  function  of  geographic  position 

Rg  Set  value  of  Earth  radius  used  for  Schuler- tuning 

S/     V  /        X  Sensitivity  of  component  or  component-group  indicated 

by  first  subscript,  with  input  and  output  indicated  by 
second  subscript 
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V.  Indicated  vertical 

(V.)  Indicated  vertical  at  start  of  problem 

V«  Reference  vertical 

r 

Vi^  True  vertical 

(V.)  True  vertical  at  start  of  problem 

W/gyx.  Angular  velocity  of  indicated  vertical  with  respect  to 

the  Earth 

^^EV^r  Angular  velocity  of  reference  vertical  with  respect  to 

the  Earth 
^^FVU  Angular  velocity  of  true  vertical  with  respect  to  the 

Earth 
Wjg  Angular  velocity  of  base  with  respect  to  inertial  space 

Wj„  Angular  velocity  of  Earth  with  respect  to  inertial  space 

W/jy\.  Angular  velocity  of  indicated  vertical  with  respect  to 

inertial  space 
cmds  Controlled  member  drive  system 

(U)  Uncertainty  in  quantity  immediately  following. 

(I)  Inaccuracy  in  quantity  immediately  following. 

1  Water  speed  compensation  current 

hfl  High  frequency  integrator 

Ifi  Low  frequency  integrator 

V^  Velocity  of  craft  through  the  water 

V^  Velocity  of  water  current 

(U)W/     V  Uncertainty  in  gyro  drift 

dir  Direct  channel  amplifier 

fb  Feedback  amplifier 

Cp  Earth  rate  compensation  voltage 

P 

Figure  1-3  represents  the  functional  diagram  of  the  instrumentation  of 
a  single  axis  of  the  complete  vertical  indicator  of  Fig.  A-1.     Figures  A -3 
and  A -4  are  the  mathematical  functional  diagrams  of  Fig.   1-3  for  the  un- 
damped and  damped  systems,  respectively. 

The  quantity  [(C) V]/^.  .^  is  defined  by  Equation  (A-1);  that  is,  [(C)V]/j.  ^v 
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is  the  angle  between  the  true  vertical  and  the  indicated  vertical. 
Then,  in  scalar  form  and  operator  notation, 

P[(C)V](t,i)=W,Ey)j    -W,Ev)i,  (A -20) 

where 

p=   -^    .  (A-21) 

dt 

Equation  (A -20)  corresponds  to  Equation  (A -2)  with  the  coupling  term 
(^(EP) Vi  ""  t(C)V](t^  i))    neglected. 

The  proposed  system  consists  essentially  of  an  accelerometer,  an  in- 
tegrator, a  single-degree -of -freedom  integrating  gyro  unit,  and  a  drive 
system  for  the  stable  platform  or  controlled  member.    The  accelerometer 
is  mounted  on,  and  near  the  axis  of  rotation  of,  the  controlled  member.    It 
will  therefore  sense  the  indicated  acceleration  of  the  base  with  respect  to 
inertial  space,  or  the  projection  of  the  resultant  specific  inertia  reaction 
force  on  the  indicated  horizontal  plane  (in  general  this  force  will  contain 
both  acceleration  and  gravity  effects). 

The  development  of  the  system  equations  is  shown  in  Derivation  Sum- 
mary A -I.    Equation  (A -32)  represents  the  solution  of  the  vertical  indicator, 
which  in  not  presented  here  as  an  analytical  expression  but  is  rather  pre- 
sented as  a  series  of  figures  in  Chapter  4,  and  is  further  discussed  there. 

INSTRUMENTATION  OF  A  SINGLE -AXIS  OF  THE  VERTICAL  INDICATOR 
OF  FIGURE  A-1  (DAMPED). 

Figure  A -4,  the  mathematical  block  diagram  of  a  damped  single -axis 
vertical  indicator,  is  essentially  a  modification  of  Fig.  A -3,  the  undamped 
vertical  indicator,  in  which  a  direct  channel  and  second  integration  with  a 
feedback  loop  are  added  in  the  orientational  control-signal-generating  loop 
to  provide  damping  and  stabilization  in  the  system.    The  essentials  of  the 
procedure  are  as  presented  in  Appendix  C. 

The  development  of  the  equations  for  instrumentation  of  this  system 
proceeds  in  an  analagous  manner  to  the  development  in  Part  II  of  this  ap- 
pendix and  is    presented  briefly  in  Derivation  Summary  A -II. 

Equation  (A  - J8)  represents  the  complete  mathematical  performance 
equation  of  the  damped  vertical  indicator.  The  solution  is  not  presented 
here  but  is  discussed  further  in  Chapter  3. 


61 


DERIVATION  SUMMARY  A-I 

The  performance  equation  for  the  accelerometer  tmit  is  (neglecting 
accelerometer  dynamics*): 

where     e    =  accelerometer  unit  outpitt  signal  , 

S/   w   N  =  accelerometer  unit  sensitivity  for  an  acceleration 
V.  "ttM  «  ;   input  and  a  voltage  output  , 

a^_  =  acceleration  of  the  "base  with  respect  to  inert ial  space  . 

Since  the  system  intrinsically  minimizes  [(C)V]/   .x,  this  angle  will 
be  small  eno-ugh  in  general  to  permit  the  assumptions: 

In  Reference  8  it  was  shown  that ,  at  submarine  velocities  » 

*IB  =  ^P^(EV)t  •  ^^23) 

Substituting  this  in  (A-.22)  gives 

•(an)  =  S(au)(a.e)  «(0)V](,_i)  +  RpW^^^j^)  ♦  (Uje^^^j  .  (^2k) 
The  integrator  performance  eqviation  is  given  "by 

^(int)  =  ^int)(e,i)  ?  (®(au)  *  ^^^^(int)  (in))  *  ^(int)  '   ^^^^^ 


*  To  incorporate  Schuler-tuning  characteristics,  the  system  is  designed 
to  have  a  natural  period  of  approximately  84  minutes.  Therefore,  the  dynamics 
associated  with  the  accelerometer  unit  and  the  controlled  member  drive  are 
neglected. 

♦•  In  this  scalar  eqriation  the  centripetal  acceleration  due  to  the  Earth's 
daily  rotation  is  included  in  g  rather  than  a__.  The  net  vertical  force  is 

effectively  g  for  submarine  operation. 
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where : 

i/j  i.\  =  current  output  of  integrator  t 

Unt; 

S/^  ^^/  !\  =  sensitivity  of  integrator  for  voltage  input,  ciirrent  rate 
(int)(e.i)   ^^^p^^^ 

^(int)(e.i)^"^^int)(in)  =  ^<=ertainty  in  i^^^^  due  to  equivalent  voltage 
input  ♦» 

i^  .X   =  integrator  current  at  start  of  problem  • 

The  primary  purpose  of  the  gyro  unit  is  to  provide  a  means  of  stalJi- 
lizing  the  controlled  mem'ber  against  ar'bitrary  motions  of  the  "base;  however, 
in  the  present  analysis,  the  gyro  unit  is  simply  considered  to  "be  part  of 
the  angular  velocity  generating  system.  The  input  current  produces  a  torque 
on  the  gyro  unit  which  is  "balanced  "by  the  angular  velocity  of  the  controlled 
meml)er  with  respect  to  inertial  space,  there'by  keeping  the  gyro  unit  on  null. 

where 

(U)M/  V  =  H(U)W/  V  =  tuacertainty  torque  on  gyro  gim"bal  , 

S/.  w.  j,x  =  sensitivity  of  gyro  unit  torque  generator  for  current 
^   »     input  and  torque  output  . 

iyf  =   current  compensating  for  rotation  of  earth  with  respect 

to  inertial  space  » 

H        =  angular  momentTom  of  gyro  \mit  » 

W/  .  V    =  angular  velocity  of  controlled  mem'ber  with  respect  to 
inertial  space  , 

W     =  W 
(IV)i    (cmds)  =  anguleo"  velocity  of  indicated  vertical  with  respect 

to  inertial  space* 
Let 


*  c.f .  Chapter  2  -  Uncertainties  in  Components 
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where 

S^  ^  w.  „\  =  sensitivity  of  controlled  member  drive 
^^     SM  ♦  ;    system  for  a  ciirrent  input  and  an  angu- 
lar velocity  output « 

From  the  foregoing, 

^IV)i  =  ^3V)i  *  h^   =  ^c,ndB)(i.W)(i(int)*  V  *  ^^^^ga)  '   ^^^S) 

i^.  =  r +   (I)i^  =  compensating  current     +     inaccuracy 

(cmd8)(i,W)  in  compensating  current  t  (A^'->q) 

\m)i  '  ^cnKl,)(i.W)('(lnt)  ♦  ">^'  *  ("'-"(gu)   •  ^^3°' 

SuTjstituting  the  above  relationships   into   Si.    (A-30)  ,  we  obtain 
^(EV)^  -  P  l^^°^^^(t,i)  =  ^(SV)t  -   ^*BV^  -  ^(EV)i^ 

=  ^cmds)(i.W)  ^*-^       [^au)(a.e)(^(°)^3(,.,)  *  V  \m)^^ 

*  (^)-M  *   (U)^int)(in)]   ^  ^cmds)(i.W)    ^^^(W)  *  ^^^^gu) 

*  ^(cmds)(i,w)   ^(int)o   *  (iU3l) 

Equation  (A-31)  corresponds  to  Eq.  (A-^).  Substituting  Eq.  (A-31)  Into 
Eq.  (A-20),  differentiating  the  result  and  setting 

^au)(a,e)  ^int)(e,i)  ^cmds)(i,W)  =  T    ^°^  ^chuler  tuning,  we 
obtain  the  differential  performance  eqiiation 
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Derivation  Summary  A-II 

The  perforaiance  equation  of  the  accelerometer  unit  is  unchanged  "by 
the  modifications  and  is: 


However,  the  output  of  the  integrator  is  now: 
g 

^mt)  '  '"""'■''   (  «(dlr)  *  «(i„t)  *  «(cp)  *   »>=(cp) 


where 


(fb)  *   ^^^®(hfi)(in))*  ^(int)o 


(dir)  -  ^(dir)(e,e)  ®(au)  *  ^^^®(dir) 


(iU2i|) 


(A-33) 


•(int)  =  ^''p^^''"^  (^an)  ^  (">^lfi)  (in))"^  «(int)c 
ipenec 

^fl))  =^fb)(i,e)  Hnt^  (^^«^ 


6/  V  =  signal  volteige  compensating  for  ship's  speed  through  the  water. 


'(fb) 


■^ — ^^^*  -^   =  performance  function  of  high  frequency  integrator t 
-^ -^  *'  ■'-      =  performance  function  of  low  frequency  integrator  . 


The  torque  "balance  equation  for  the  gyro  xmit  remains  uncharged, 
^^^^(gu)  *  S(tg)(i.M)  (^int  *  V  =  ^^cmds)  =  ^(IV)i  ' 


Letting  S 


(cmdsXi.W) 


'(tfi)(i.M), 


we  obtain 


^(IV)i  -  ^(E7)l  *  ^(IB)  =  ^^(Icmds) 


(A-26) 
(A-27) 
(A-28) 
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i^=S *   (Di^.  •(A-29) 

(cmde) (i,w) 


Substituting  the  above  relationships  into  Eq.  (A-30)  gives 
(P  -^  ^hfi)(e,i)  ^fb)(i.e)^  <^EV)^  -  P  C^°)^J(t,i)^ 

=  ^(hfi)(e,i)  ^(cmd8)(i.W)(  ®(dir)  *  ®(int)  *  ^(cp)  "*"  ^^^®(cp) 
-  (U)e(^j  *  We(^,^j^^^jy((U)W(^j  .  (I)i^  ^cMs)(iW)) 
KP*^hfi)(e,i)   ^fb)(i.e)) 

=  S(hfi)(e.i)  S(cMs)(i.W)[(5(°)^J(t.i)  *^V)J 

xrs  s  +     ^(lfi)(e.;)   ^(au)(a.e)N 

''l^dir)(e,e)  ^au)(a,e)  *  p  ) 

V   (dir)(e,e)  p  /  "   '    (au)  p  (lfi)(in) 

^  (">«(dir)  *  ^cp)  *  (^^^cp)  -  ^^^^fb)  ^  ^U>^hfi)(in)*  ^int)o] 

<P  *  ^hfi)(e.i)  ^fb)    (i.e)X^)^gu)  *  ^cMs)(i.W)(^)^)- 

(A.3^) 

^®**^"^  ^(hfi)(e.i)  ^(cinds)(i,W)  ^(au)(a,e)  ^(dir)(e,e)  "  r"  ^°'' 

'         (A-35) 
Schuler  tuning  and  differentiating   (A-34)  and  after  re-arreoiging,  we  get 

=  [p-(l  -  5->  P  (S(,,,j(^^,j   ^fb)(i.e)  -  H^S^^^^^^^^^^^JNeV), 
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(U)e, 


_l(au)___/   ^  °(lfl)(e.e)\    ^       .  ^  ,  ,^. 

■  \   ^au)(a.e)  ^'   ^dir)(e.e)^"  "^  ^hfi)(e.i)  ^cn>ds)  (i.w)(  ^^^^dir) 


^(hfiXe.i)  ^(cind8)(i,w)  ^(lfi)(e,e)^^^®(lfi)(in) 


-  (P^  *  P  ^hfi)(e.i)  ^n))(i,e)j  (^'^^(gv)  -  ^cmd«)(i.w) 
With  R  =  R  very  nearly,  Wrigley    has  shown  that  the  temfl  -  g-jwill  "be 


negligibly  small  under  conditions  which  will  "be  encountered  in  service.  In 
Appendix  (C),  it  is  shown  that  it  is  possi"ble  to  miniBiize  the  term  in 

s       * 
r[S(hfl)(e.l)  ^n.)(l.e)  -  s|^''jj;;°i]  «W^  ^"  ^  ^1'-  ^^ 

ratio. 

where  V  =  component  of  craft  velocity  with  respect  to  the  water  perpen- 
dicular to  input  axis  of  system., 

and  V  =  component  of  ciarrent  velocity  with  respect  to  the  earth  perpen- 
dicular to  input  axis  of  system. 

/^hfi)(e.i)^(ft)(i.e)  -  S...    ..       .\V^ 

Now,setting  e.     .  =     — g- —^ (dirji^^  \w  ^  (^3^) 

^cp;   y   S(hfi)(e,i)  ^(cmd8)(i,w)        r- 


the  terms  in  e/  \   and  W/__v  disappear  from  Hi.  (A-36)  leaving  only  a  term 
t 
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will  In  general  "be  negligitle,  "because  V  ,  the  rate  of  change  of  current 
velocity,  is  in  itself  generally  small,  and,  in  addition,  because  the  coef- 
ficient of  V  is  optimized  (Chapter  4),  Note  that  the  effect  of  any  errors 
in  the  measurement  of  (V  )  are  considered  In  the  ahove  system  equations  in 
the  term  (I)e/  v. 

Equation  (A-36)  represents  the  differential  performance  equation  for  the 
system.  .If  now  the  Laplace  transform  of  the  performance  equation  is  taken, 
the  initial  conditions  appear  as  follows  where  now  p-»-Laplace  operator. 

-r.  fs  •  s  ^(lfi)(e.;)/c 

-n^hfi)(e.i)^fb)(i.e)-S(^^^j(^^^);R 

^    ^(dir)(e,e)  ^  ^s  ^(au)(a,e) 
-  P  ^hfi)(e.i)  ^cxads)(i.w)  [^"^^dir)  *  ^^^^cp)  "  ^^^^ft) 

*  ^^^®(hfi)(in)]  -^(hfi)(e,i)  ^(cmd8)(i,w)  ^(Ifi)  (e.^)  ^^^®(lfi)  (in) 

-(P^  *PS(hfi)(e,i)  ^n>)(i.e))  ^^^V)  -  ^cmds)(i.w) 
^  (P^  *  P  ^hfi){e.i)  ^ft)(i.e))  ^'K  *  C(C)Vj(,,i), 
H^*'(hfi)(e.i)  ^fl.)(i.e))  ^^'^''ht,i)o^{^'^^^(Mi)(e,i)   ^(fb)(i,e) 

*i:)^^°)^J(t.i)o-(^hfi)(e,i)3cf.)(e.i)-^^)'H^ 

(n)e/  V 

*  R  S^auK^.e)  *  ^(hfi)(e.i)  ^(cmds)(i,w)((^^^dir)o  *  ^^^^cp)o 
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"(^  *  ^(hfi)(e,i)  ^(fb)(l.«))    ^"^^(gu)o  *  ^(cmds)(i,w)    ^^^^(w)o 

*  ^(cmds)(i.w)  (P  *  ^hfi)(9.1)   ^n,)(i.e))^^^^(w)o'  (A-38) 
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APPENDIX  B 

POSITION  INDICATION  BY  OPEN-CHAIN  INTEGRATION 
OF  THE  INDICATED  VERTICAL 


Figures  B-1  and  B-2  are  the  mathematical  functional  diagrams  of 
the  proposed  position  indicators  corresponding  to  the  damped  and  un- 
damped vertical  indicators  respectively,  as  described  in  Appendix  A. 

Derivation  Summary  B-I  presents  the  mathematical  derivation  of 
position  indication  as  derived  for  both  systems.    This  derivation  has 
the  same  form  whether  the  vertical  indicator  is  damped  or  undamped. 
This  fact  is  readily  apparent  when  it  is  recalled  that  the  relationship 
between  the  angular  velocity  of  the  controlled  member  and  the  inte- 
grator current  is  the  same  regardless  of  damping.    The  relationship 
is  given  by  Equation  (A  -30) . 

Equation  (B-6)  shows  that  the  correction  to  the  indicated  position 
is  the  sum  of  the  correction  to  the  indicated  vertical,  initial  conditions 
in  position  and  correction  to  the  vertical,  together  with  the  time  integral 
of  certain  uncertainties  and  inaccuracies. 

Recalling  that  the  transfer  function  has  been  defined  as  the  response 

of  the  system  to  the  forcing  function  e^^f  ,  we  may  write    F(w,  t)|  for 

the  open-chain  integrator  as   -— —  (1  -  e-'^f  )  .    The  mean-square  er- 

I    ^f  I 

ror  may  then  be  evaluated  quite  dimply  after  the  method  applied  to  the 

undamped  closed-chain  vertical  indicator  in  Chapter  3.    That  is,  the 

plot  of  1-.^^^  (1  -  e-l^f  )  I  2  may  be  constructed  and  the  integral  of  Equation 

(2-11)  apiAied. 


^(output) 


-.=/' 


F(w,t) 


2     G(.^^(w)    dw,  (2-11) 


where  G^^^^  represents  the  power  spectral  density  of  each  of  the  uncer- 
tainties shown  in  Equation  (B-5)  taken  in  turn. 
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Derivation  Svunmary  B-I 


Position  Indication  "by  Open- Chain  Single  Integration 
of  the  Indicated  Vertical 


In  both  the  systems  of  Figs.  B-1  and  B-2,  position  information  is  ob- 
tained from  a  single  integration  of  the  current  input  to  the  gyro  unit. 
Writing  the  differential  performance  equation  for  the  open  chain  integration: 

where   ^^^^^piWin^  ~  i"icertainty  in  output  of  the  position  indicating 
integrator  reduced  to  an  equivalent  signal  input, 

(P.)  =  initial  position. 

1  0 

FromEq.  (A-30),and  noting  that  W^^j^^^  =  W^^j^^^  -  p[(C)V]^^  ^^  from  Eq.  (A-.20), 

^(-)=  %;z:^  (V)t-pc(c)v],,,,,-(u).v,)-(i)i.. 

(B.2) 

Substituting  Bq.  (B-2)  into  Bq.  (B-l)  and  arbitrarily  setting 

^(cM.)(i.w)=^Pl)(i?)'  '"°'""'' 

-^(PiXi.P)  "''w)*<^'o  •  (=-»> 

Noting  that  I  W(j^j^  =  Pj  -(Pj)^  .  we  write 
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*  ?((«^(Pl)(ln)  -   «"«(^)  -  ^H)(l.P)   "H,  )   .  <"> 

IHirther .defining    [{0)Pj,|.  i)  =  I'l  -  ^i  . 

ccWd.i)  =  i:(°W(t,i)o  *  c«=)^J(t.i)  -  [(=>^:i(t.i). 

Writing  this  in  the  more  conventional  form  we  get 

-  •fo((''>^Pi)(in)  -  »>''(«u)  -  ^Pi)(lP)  "Hj"-  <=-«) 
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APPENDIX  C 


OPTIMIZATION  OF  SYSTEM  PARAMETERS 


In  any  dynamical  system,  the  engineer  is  forced  to  a  compromise  be- 
tween a  system  with  minimum  forced  or  steady- state  errors  and  a  satis- 
factory transient  response. 

Given  a  performance  function  of  the  form 

Blasingame^  '  has  concluded  that  an  optimum  practical  system  may  be 
achieved  with  a  performance  function  whose  characteristic  equation  is 
a  cubic;  and  furthermore,  such  a  system  may  be  achieved  by  a  network 
involving    only  integrators  and  summing  devices,  as  shown  in  Fig.  C-1, 
where  the  C's  and  S's  may  have  all  values,  positive  or  negative,  including 
zero. 

Such  an  optimum  system  has  been  applied  to  the  undamped  vertical 
indicator  of  Fig.  A- 3  to  obtain  the  damped  vertical  indicating  system  of 
Fig.  A-4. 
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Fig.  C-1    Generalized  mechanization  computer. 
The  performance  equation  for  the  vertical  indicator  of  Fig.  A- 4  for 
an  input  in  the  angular  velocity  of  the  true  vertical  may  be  written 


[(C)V],,)        ^^'"V)(a,W,H)-P(^-^^-^ 


w 


(EV)t 


P^2(DR)W^3P.W2^3    p.(^ 


(C-2) 


P^(l-S(Vi)(a,W)^-P^d 


,A_  .S(Vi)(a„W)^ 


•"^^t^^^t-tt 


where 


(C-3) 


^(Vi)(a,W) 


S^"(a,e)^(dir)(e,e)^(hfi)(e,i)^(cmds)(i,W) 


1 
-^d 


ns      (Cl)j)         ^(dir)(e,e) 


(C-4) 
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W2 

ns 


2      (CT)      =    S(i^i)(g  i)S(j^)(.  g)    ,  (C-5) 

"  nE  ^ 

2(DR)W„^ 
(CT)j5  =  -^ f"    '  (C-6) 

W^g    =   natural  angular  frequency  of  the  system, 

^nE    =   earth's  natural  angular  frequency, 
(CT)j^   =   system  characteristic  time. 

2 
The  p    term  in  the  numerator  of  Equation  (C-2)  has  the  significance 

of  an  acceleration  in  the  angular  velocity  of  the  true  vertical  and  gives 
rise  to  errors  known  as  jerk-errors  in  the  system.  The  p  term  in  the 
numerator  of  Equation  (C-2)  has  the  significance  of  a  rate  of  change  in 
the  angular  velocity  of  the  true  vertical  and  gives  rise  to  errors  in  the 
system  known  as  velocity  errors. 

The  physical  significance  of  the  two  changes  associated  with  the  angu- 
lar velocity  of  the  true  vertical  may  be  seen  if  it  is  recalled  that  movement 
of  a  craft  over  the  surface  of  the  earth  produces  a  change  in  the  direction 
of  the  true  vertical.    This  change  in  direction  or  angular  velocity  of  the 
true  vertical  may  be  represented  as  a  vector  W/^yv    at  right  angles  to 
the  radius  vector  of  the  earth  and  the  velocity  vector  of  the  craft,  and  of 


magnitude 


ZSL 


Re 


,  where 


is  the 


is  the  velocity  of  the  craft  and  Rg 
radius  of  the  earth. 

Thus,  a  constant  linear  acceleration  of  the  craft  over  the  earth  re- 
sults in  a  steady  rate  of  change  in  the  magnitude  of  the  vector  W/^y^  , 
and  is  physically  associated  with  the  p  term  in  the  numerator  of  Equa- 
tion (C-2).    However,  maneuvering  accelerations  normally  encoimtered 
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on  shipboard,  say  in  a  change  of  speed  from  5  to  15  knots,  involve  a  rapid 
change  of  acceleration  from  zero  to  a  maximum  of  about  5  knots  per  min- 
ute, with  another  rapid  drop  in  acceleration  as  the  new  speed  is  approached. 

2 
These  changes  in  the  rate  of  acceleration  are  associated  with  the  p    term 

in  the  numerator  of  Equation  (C-2),  and  numerically  are  considerably 
greater  than  the  maximum  accelerations  involved.   Similarly,  a  change  of 
course  involves  a  change  in  the  magnitude  of  the  components  of  W/^yv 
along  each  of  the  two  horizontal  axes  of  the  vertical  indicator  that  is  sta- 
bilized in  azimuth.   Thus,  a  course  change  also  gives  rise  to  both  jerk 
and  acceleration  errors,  in  which  the  rate  of  change  of  acceleration  will 
again  exceed  the  acceleration  for  a  normal  tactical  maneuver. 

An  examination  of  Equation  (C-2)  shows  that  it  should  be  possible  to 
eliminate  one  of  the  terms  in  the  numerator  by  a  proper  choice  of  sensi- 
tivities, and  then  to  minimize  the  effects  of  the  remaining  term.  A  study 
of  acceleration  and  jerk  effects  likely  to  be  encountered  on  shipboard  in- 
dicates that  it  may  be  desirable  to  proceed  on  the  basis  of  eliminating  the 
jerk-error  term  and  minimizing  the  acceleration- error  term,  as  has  been 
done  in  the  following  sections. 

With  the  jerk-error  term  eliminated  by  setting  S/y.w    w)   =  p  > 
Equation  (C-2)  may  be  written 


Letting 


^(EV),    -     p2,2(DR)W,3.w2^3(p.^) 


pjWf   =   j(FR)W„E    , 

Wj  =  angular  frequency  of  forcing  function , 
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(C-7) 


(FR)   =    ^    , 


and  substituting  into  Equation  (C-7),  we  obtain 

,2 


3(FR)W„^ 


4(DR)^(FR)^^+(l-(FR)^^^ 


1  -  (FR)2^ 


(EV)t  W  nE^^^^  n    /     (FR)^  ^ .  ^ .  2(PR)(rR)V.  .frn^     (ct^    + 1 

(C-8) 
Hydinger  and  Webster ^^  have  shown  that  the  forcing  frequencies  in 

W/gyx     likely  to  be  encountered  on  shipboard  will  be  essentially  of  zero 

frequency  or  such  that  (FR)  >  10. 

Optimizing  on  the  basis  of  (FR)  >  10,  it  may  be  shown  that  the  following 

approximations  hold: 


a.  (FR)W^e(CT)j)  >>1 


^■^)»^^^^>^ 


By  making  use  of  these  approximations,  Equation  (C-8)  is  reduced  to 
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t(C)V](t.i) 


4(DR)^FR)^„-f(l-(FR)^„y 


1  -  (FR)' 


(FR)W„E 


W 


(EV). 


<E^^<  [-  ^j(i  (FR)W„e(CT)d 


(C-9) 


(ct)j5w2j^e(fr)2 


2 
4(DR)^(FR)^^-^(l-(FR)\)   \ 


l-(FR)' 


(C-10) 


Substituting 
(CT)jj 


2(DR)W, 


ns_ 


2(DR)(FR), 


WnE-<s    '   W„Etl-(FR)'n] 


(C-6) 


and  rearranging,  we  obtain 


[(C)  V] 


w 


iUl  _ 


X 


(EV). 


2(FR)^W, 


nE 


^^""^^^^^^n-imm. 


2(FR)^-(FR); 
(DR) 


(C-11) 


Assuming  a  fixed  (DR)  and  taking  the  partial  derivative  of  Equation 
(C-11)  with  respect  to  (FR)_,  we  have  the  necessary  condition  for  a 


minimum  in 


[(C)V] 


(U) 


W/ 


(EV)^ 


to  exist  at  a  given  damping  ratio. 


W 


(EV) 


2(FR)^W„^ 


4(FR)^^(DR)^  -  1  -  2(FR)^^-h  3(FR)^ 
(FR)2(DR) 


(C-12) 
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Solving  Equation  (C-12)  for  (DR)  as  a  function  of  (FR)^^,  we  have 

2         1  +  2(FR)^.3(FR)^ 

(DR)2    =  B— ^ a-       .  (c-13) 

4(FR)2 

Figure  C-2  is  a  plot  of  (DR)  vs.  (FR)^^  from  Equation  (C-13).    With 

(FR)j^   =    1,  (DR)    =   0,  and  with  these  conditions  in  Equation  (C-13),  we 

[(C)V],^.. 
find  -TTj ^-*-^   =  0.    The  conditions  are  essentially  those  of  an  undamped 

<^A  (13) 

Schuler-tuned  inertial  system,  which  Wrigley^     '  has  shown  to  be  insensi- 
tive to  a  sinusoidal  forcing  function  with  an  84. 4- minute  period  correspond- 
ing to  (FR)jj   =    1. 

However,  a  condition  of  zero  damping  and  (FR)      =    1  is  neither 
attainable  nor  desirable  in  practice  for  an  inertial  system  operating  over 
an  extended  period.    Hence,  choosing  (DR)  =  0.1,  0.3,  and  0.5  for  a  com- 
parative study,  and  utilizing  the  relationships  of  Equations  (C-16)  through 
(C-18)  as  shown  below,  we  arrive  at  Table  C-1,  which  is  a  tabulation  of 
the  system  parameters  for  the  vertical  indicator  of  Fig.  A-4,  as  obtained 
when  optimum  values  of  (FR)    are  chosen  for  the  various  damping  ratios. 
The  resulting  characteristic  equations  are  also  shown  below. 

The  characteristic  equation  of  (C-2),  which  is  in  the  form  of  a  cubic, 
can  also  be  expressed  as  the  product  of  a  second-order  term  and  a  first- 
order  term  as  follows: 

(p^  +  agP^  +  a^  p  +  a^)    =    (p^  +  Ap  +  B)(p  +  D)  ; 

(C-14) 
also  as 

(p2  +  AP  +  B)(p+D)    =    (p2  +  2(DR)W^gP  +  w2g)(p  +  ^^)    . 

(C-15) 
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Equating  coefficients  of  like  terms,  we  get 


ag   =    (A  +  D)    =   2(DR)W^g+^^^j       ,  (C-16) 


9        2(DR)W„, 
^1    =    (B  +  AD)    =   W^g  +       (CT)"      '  (C-17) 


w2 


2(DR)W  2(DR)(FR)_. 

Figures  C-3  and  C-4  are  the  resulting  steady -state  responses  of 
the  optimized  vertical  indicator  to  sinusoidal  forcing  motions  of  the 
true  vertical.   It  will  be  noted  that  the  response  to  accelerations  is 
three  orders  of  magnitude  greater  than  the  response  to  angular  velo- 
cities at  (FR)  =  1;  however,  the  response  to  accelerations  is  decreas- 
ing at  18  decibels  per  octave,  whereas  the  response  to  velocities  is  de- 
creasing at  12  decibels  per  octave.   Consequently,  at  (FR)  =  100,  the 
responses  to  both  inputs  are  of  equal  magnitude,  while  at  higher  fre- 
quencies the  response  to  accelerations  becomes  increasingly  lower 
than  the  response  to  velocity  inputs. 
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APPENDIX  D 


METHOD  USED  IN  THE  DETERMINATION  OF  SYSTEM 
PERFORMANCE  CHARACTERISTICS 


Consider  a  generalized  performance  equation  of  the  form 

(p3  +  agP^  +  a^p  +  a^)  q(p)^^^   =   q(P)(in^)H-Pq(P)(in2)  +  (P  +  ao)*l(P)(in3) 

^-P(p  +  ai)q(p)(i^^)  -  <l(in2)<°+)  '  ^(ing)^^^)  '  <P  +  ^l^^Cin^)^^) 

-  ^(in4)<^) »  CD-I) 

where 

p  =  Laplace  transform  variable   y 

q(p)^^  =  Laplace  transform  of  q(t)jjj    , 

q/jjjx(o+)  =  initial  value  o^  Q/jn)  at  t  =  (o+),  and 

a  ,  a^,  a^.  .  .  are  numerical  constants. 

The  complete  response  of  the  differential  equation  may  be  found  as  the 
sum  of  the  responses  to  the  various  terms  on  the  right-hand  side  of  the 
equation.   Thus,  the  generalized  performance  equation  may  be  broken  down 
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into  a  series  of  performance  equations  with  but  a  single  term  on  the  right, 
such  that  q(p)Qyt.  may  be  represented  as  follows: 
Setting 

[P]    =  (  p    +  agP    +  a^p  +  a^)  ,  the  characteristic  equation, 
we  obtain 

^(P)(outi)    =  Tp]^^)(inp    '  (D-2) 

<l(P)(out2)    =  rF]q(P)(in2)    '  (D-3) 

^(P)(out3)    =    -TFT  ^^)(in3)    '  <D-4) 

q(p)(out  )  = and  so  forth. 

Expressions  (D-2),  (D-3),  and  (D-4)  show  that  the  generalized  system 
may  be  considered  as  made  up  of  a  number  of  independent  sub-systems 
with  performance  equations  given  by 

(p  +  aj 


JL 


O' 


-«  2  '32  '  ^         3  2 

p'^  +  agP^  +  a^p  +  a^  p  +a2P  +a^p  +  aQ  p  n-agP  +2i^p+SL^ 


respectively.      Thus,  the  weighting  fimction  for  the  generalized  system 

for  an  input  in  q(t)/^  \  may  be  determined  quite  simply  by  setting  the  sub- 

^    n' 
system  of  q(t)/.j^  .  into  an  analog  computer  and  taking  the  time  response  to 

^    n' 
a  imit  impulse. 

It  should  be  noted  here  that  the  terms  which  have  to  do  with  the 

initial  conditions  are  in  the  form  of  a  product  of  a  performance  equation 

and  a  constant,  and  that  the  inverse  Laplace  transform  of  a  constant  has 

the  form  of  an  impulse;  hence  the  time  response  to  initial  conditions  has 
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the  form  of  a  weighting  function.   Similarly,  if  qW/^jj)  has  some  probable 
form  that  may  be  predicted  on  the  basis  of  its  physical  nature,  the  probable 
response  of  the  physical  system  to  q(t)Qj.\  may  be  found.   In  particular, 
if  the  system  inputs  and  outputs  are  represented  as  in  the  mathematical 
block  diagram  of  (D-5), 


(D-5) 


JWft 

afF(Wf,t)e 

¥ 

a^e 

(PF)  system 

then 


where 


jw^t  I  I    J(Wft  +  A,     ) 

^^t)(out)    =   afF(Wf,t)e     ^     =   aj|F(Wf,t)|e      ^        ^y" 


(PF)g      =   system  performance  equation  of  form  of  (D-2)  etc., 

F(Wf,t)   =  transienttransf  erf  unction  f  or  frequencyw.  (c.f. 
^  Chapter  2),  ^ 

A  .  =  dynamic  phase  angle  associated  with  F(w,t)    ; 

and,  if  it  is  desired  to  find  the  absolute  value  of  I(w£,t)  as  a  function  ci 
time  for  a  given  frequency  w-,  then  the  output  of  the  system  to  the  forcing 

function 

jw,t 
a- e  =   a£(cos  Wxt+jsin  Wjt )   ,  (D-6) 

for  which  the  system  response  may  be  written 

^^*)(out)    =   aj|F(Wf,t)|cos  (Wjt+A^y^)  +  jaf|F(Wf,t)|  sin(Wjt  +  A^yj^)    , 

(D-7) 
may  be  represented  on  the  complex  plane  as  in  Fig.  D-1. 
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Axis 
Imagina 


Axis  of 
Reals 


a^ |F(wpti)|cos  (Wft^  +  A^y^) 

Fig.  D-1 
Complex  plane  representation  of  F(wpt)e^^f  . 

From  Fig.  (D-1),  it  is  apparent  that 


Solving  for    F(wj,ti  we  get 

I  ^(*)out 


F(Wf,t)| 


(D-8) 


(D-9) 


1  **f  *■  "  "l 

The  mathematical  relationships  of  the  foregoing  paragraphs  suggest 
that  |F(wj,t)|  may  be  obtained  quite  simply  by  use  of  a  high  speed  analog 
computer  utilizing  timing  marks  to  indicate  time  and  clamping  to  ensure 
zero  initial  conditions  at  the  start  of  each  computing  cycle. 

If  the  input  to  the  system  is  chosen  as  a  sinusoidal  input  of  unit 
amplitude  such  as  cos  w^t,  the  output  will  be  obtained  as  a  function  of 
time  such  that 


q(t) 


out  =|F(Wf,t)|(^°^V  +Vn 


(D-10) 


If  then  the  output  (D-10)  is  put  through  an  integrator  with  the  performance 
function  w^/p ,  the  result  will  be 


UHCLASSIFIED 


q(t)ou^  =    F(wpt)w^  J    coswjt 
0 


dt 


=  I  F(wpt)  I  sin  (wjt  +  A^y^).  (D-11) 


If  the  output  (D-10)  is  now  taken  to  the  x-axis  input  of  an  oscilloscope 
and  the  output  (D-11)  is  taken  to  the  y-axis  input  of  the  oscilloscope, 
the  picture  on  the  oscilloscope  will  be  the  time  varying  presentation 
of  Fig.  (D-1),   where  the  radius  from  the  center  of  the  diagram  to  the 
resulting  curve  is  representative  of  the  sum  in  Equation  (D-7),  which 
has  been  shown  in  Equation  (D  -9)  to  be  equivalent  to 


F(wpt)| 


q(t)out 


or 


F(wpt)      =   q(t)Q„4.      =  radius  of  oscilloscope  curve  at 

time  (t),  for  a  sinusoidal  input  amplitude. 

The  function  |  F(w,t)|,which  was  derived  in  Chapter  2,  may  thus  be 
determined  by  applying  unit  forcing  functions  of  various  frequencies 
over  a  frequency  band  for  which  the  response  is  desired  and  taking  the 
amplitude  of  the  radius  of  the  oscilloscope  presentation  as  a  function  of 
time. 

It  is  apparent  from  Fig.  D-1  that  knowing  the  frequency  and  the  time 
for  any  point  on  the  curve,  and  having  a  base  line  established  by  some 
means,  the  dynamic  phase  angle  may  be  easily  determined  from  the 
relationship 

A(dyn)    =A-Wft. 
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